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ABSTRACT 


This  report  describes  the  heating  properties  of  the  Canadian 

Forces  Armoured  Personnel  Carrier  (APC).  In  the  theoretical  analysis,  the 

APC  Is  modelled  as  a  box,  having  various  Insulating  covers  on  the  surfaces. 

The  rate  of  heat  transfer  from  the  APC  Is  then  estimated  for  each 
Insulation  configuration  using  empirical  heat  transfer  coefficients  and 
measured  thermal  conductances.  All  assumptions,  equations  and  physical 

data  used  In  the  analysis  are  presented  as  Is  a  listing  of  the  FORTRAN 
computer  program  which  was  used  to  solve  the  heat  transfer  equations. 
Thus,  the  reader  may  use  this  report  as  a  guide  to  performing  the  heat 
transfer  analyses  for  similar  problems.  Field  measurements  of  Internal  air 
temperature  and  of  heat  flow-rates  were  made  using  two  APC's  In  various 
Insulating  configurations.  The  measured  heat  loss  from  an  uninsulated  APC 
was  In  close  agreement  to  the  heat  loss  predicted  by  the  theoretical 
analysis.  However,  the  measured  heat  loss  from  an  APC  Insulated  with  a 
covering  developed  by  the  Defence  Research  Establishment  Suf field  was  found 
to  be  significantly  greater  than  that  predicted  by  the  theoretical 
analysis.  Theories  are  advanced  on  possible  sources  of  error  and 
discrepancy  between  the  model  predictions  and  the  measured  values  of 
temperature  and  heat  flow.  Heat  loss  due  to  ventilation  Is  also  examined 
briefly. 
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RESUME 


II  s'aglt  d'un  rapport  decrlvant  les  proprietes  thermlques  du 
vehlcule  bllnde  de  transport  de  troupes  (VBTP)  des  Forces  canadlennes. 
Dans  1 'analyse  theorlque,  le  VBTP  est  represente  par  une  bofte  recouverte 
de  dl verses  enveloppes  isolantes.  On  estlme  ensulte  le  taux  de  transfert 
de  la  chaleur  qul  se  degage  du  vehlcule  dans  chaque  cas,  utlllsant  i  cet 
effet  des  coefficients  emplrlques  de  transfert  de  la  chaleur  et  les 
conductances  thermlques  mesurees.  Toutes  les  hypotheses,  les  equations  et 
les  donnees  sur  les  proprietes  physiques  utlllsees  aux  fins  de  1'analyse 
alnsl  que  le  programme  d'ordlnateur  en  langage  FORTRAN  ayant  servl  a 
resoudre  les  equations  de  transfert  thermlque  y  sont  presentes.  Alnsl,  ce 
rapport  peut  servl r  de  guide  au  lecteur  pour  effectuer  des  analyses  de 
transfert  thermlque  appllcables  3  des  problemes  slmllalres.  La  temperature 
de  1‘alr  a  l'lnterleur  de  deux  VBTP  recouverts  de  dlverses  enveloppes 
Isolantes  alnsl  que  le  flux  de  chaleur  qul  se  degagealt  des  v§h1cu1es  ont 
ete  mesures  en  situation  reelle.  La  perte  de  chaleur  du  VBTP  non  Isole  tel 
que  mesuree  Concorde  tres  blen  avec  les  previsions  theorfques.  Par  contre, 
11  y  ent  un  ecart  considerable  entre  les  resultats  experlmentaux  et  le 
modele  theorlque  dans  le  cas  du  VBTP  avec  l'enveloppe  Isolante  du  centre  de 
recherches  pour  la  defense,  Suffleld.  Des  theories  sont  avancees  sur  les 
causes  possibles  d'erreur  et  les  divergences  qul  existent  entre  les 
predictions  theorfques  et  les  mesures  obtenues.  La  perte  de  chaleur 
attrlbuable  a  la  ventilation  est  egalement  examinee  en  bref. 
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LIST  OF  VARIABLES 


A,A.«  -  surface  area. 

[A]  •  -  a  matrix  whose  elements  are  composed  of  the  various  thermal 

resistances  of  the  enclosure. 

A(1,j)  -  elements  of  the  matrix  [A]. 

[B]  -  a  vector  whose  elements  are  composed  of  the  known  air 

temperatures  and  the  thermal  resistances  of  the  enclosure. 

B(i)  -  an  element  of  the  vector  [B]. 

c  -  an  empirical  heat  transfer  constant. 

Cq  -  fluid  heat  capacity  at  constant  pressure. 

f!|  j  -  the  viewf actor  from  surface  1  to  surface  j. 

g  J  -  the  acceleration  due  to  gravity. 

Gr  -  the  Grashof  Number,  gpL  ^Tj-TJ/v2. 

h  -  the  convective  and  conductive  heat  transfer  coefficient. 

1,j  -  Integer  variables  designating  surfaces, 

k  -  fluid  thermal  conductivity. 

L  -  the  characteristic  length  of  a  surface, 

n  -  an  empirical  heat  transfer  constant. 

Nu  -  the  Nusselt  Number,  hL/k. 

Pr  -  the  Prandtl  Number,  Cpv/k. 

q  -  the  fluid  volumetric  flowrate. 

Q  -  the  total  rate  of  heat  loss  from  an  enclosure  from  all  sources. 

Qj  -  the  net  radiant  heat  transfer  rate  leaving  surface  1. 

Qv  -  the  rate  of  heat  loss  from  an  enclosure  due  to  air  exchange  with 
the  environment. 

r^j  -  the  radiative  thermal  resistance  between  surface  1  and  surface  j. 

Re*  -  the  local  Reynold's  Number  based  on  position,  Ux/v. 

Re,  -  the  Reynold's  Number  based  on  the  characteristic  length,  UL/v. 

Rea  “  the  equivalent  thermal  resistance  of  the  enclosure. 

R^j  -  the  jth  convective  or  conductive  thermal  resistance  of  surface  1. 

R'.j j  -  the  net  jth  thermal  resistance  of  surface  1  due  to  the  action  of 

J  convective  and  radiative  thermal  resistance  In  parallel. 


LIST  OF  VARIABLES  (Cont'd) 


St  -  Stanton  Number,  h/(pC_U). 

t  -  thickness.  M 

U  -  measured  local  Internal  air  temperature. 

[T]  -  a  vector  whose  elements  are  the  desired  internal  wall  surface 

temperatures. 

Ta  -  the  volume  weighted  internal  air  temperature. 

T.j  -  an  element  of  the  vector  [T]. 

Tm  -  a  mean  temperature. 

T0  -  the  internal  vehicle  air  temperature. 

T„  -  the  external  ambient  air  temperature. 

U  -  the  free-stream  fluid  velocity  over  a  surface. 

V  -  the  total  volume  of  an  enclosure. 

V.j  -a  local  volume  element  associated  with  the  temperature  tj. 
x  -  the  length  from  the  leading  edge  of  a  flat  plate  to  the  point  of 

interest. 

p  -  fluid  coefficient  of  thermal  expansion, 
e  -  the  emlssivity  of  a  surface, 

v  -  the  fluid  kinematic  viscosity, 

p  -  the  fluid  density. 

a  -  the  Boltzman  Constant,  5.669  x  10"8  [W/m2K4]. 


( vl  1 ) 


1.0  INTRODUCTION 


The  Canadian  Forces  (CF)  have  a  commitment  to  operate  In  regions 
In  which  the  weather  can  be  as  much  a  problem  as  any  human  adversary.  The 
CF  are  frequently  on  exercise  In  Canada's  High  Arctic  and  in  northern 
Norway  where  temperatures  can  easily  reach  -40#C.  Problems  associated  with 
the  environment  can  often  be  overcome  by  making  use  of  specialized 
equipment  and  relevant  training. 

Mechanized  Infantry  units  of  the  CF  are  using  Armoured  Personnel 
Carriers  (APC's)  to  Increase  their  mobility  and  effectiveness. 
Unfortunately,  under  the  severe  cold  weather  conditions  which  can  occur  In 
the  Arctic,  vehicle  problems  are  abundant.  In  an  attempt  to  reduce  these 
vehicle  problems,  conserve  fuel  and  keep  the  occupants  of  the  APC's  warm. 
Defence  Research  Establishment  Suffleld  (DRES )  was  tasked  to  devise  an 
Insulating  cover  for  the  APC.  The  primary  goal  of  the  task  Is  to  produce  a 
cover  which  will  maintain  the  Internal  vehicle  temperature  at  a  level  that 
will  facilitate  starting  of  the  vehicle  In  extremely  cold  weather,  with  a 
minimum  of  auxiliary  heating. 

DRES  approached  the  Environmental  Protection  Section  (EPS)  of 
Defence  Research  Establishment  Ottawa  (DREO)  for  advice  on  the  thermal 
properties  for  a  protective  covering.  The  EPS  Is  devoted  to  research  In 
the  field  of  protecting  the  man  from  various  aspects  of  the  environment, 
but  specializing  In  cold  weather  protection.  Protecting  a  man  from  the 
cold  has  many  similarities  to  the  problem  at  hand. 

This  report  presents  an  Idealized  analysis  of  an  APC  which 
provides  the  heat  transfer  Information  which  DRES  requires  to  choose  Its 
most  effective  course  of  action  In  developing  a  vehicle  covering.  This 
report  gives  the  reasoning  and  procedures  behind  each  step  or  assumption  of 
the  analysis.  Thus,  the  reader  may  use  this  report  as  a  guide  for 
performing  similar  analyses  of  heat  transfer  from  enclosures.  With  a 
minimal  amount  of  work,  the  Interested  reader  can  perform  the  analysis 
under  different  sets  of  ambient  conditions,  with  different  Insulation 
values  or  under  different  sets  of  assumptions  as  required  by  particular 
problems.  The  algorithm  used  to  perform  the  actual  calculation  of  the  heat 
transfer  and  temperatures  (Appendix  C)  Is  written  In  standard  FORTRAN-77 
and  was  used  on  a  Honeywell  computer  under  a  CP-6  operating  system.  For 
the  layman,  only  a  basic  understanding  of  FORTRAN  and  computers  Is  required 
In  order  to  apply  the  heat  transfer  analysis  to  similar  applications.  The 
equations  used  to  determine  the  heat  transfer  coefficients  are  given 
through  the  text  In  their  appropriate  sections. 


Experimental  measurements  of  the  effectiveness  of  the  DRES 
prototype  covering,  were  made  at  DREO  and  are  Included  In  this  report.  The 
purpose  of  these  measurements  Is  mainly  to  validate  the  theoretical 
estimate  of  the  predicted  effectiveness  of  an  insulated  cover.  A  secondary 
function  of  the  experimental  measurements  Is  to  obtain  a  quantitative 
evaluation  of  the  DRES  prototype  covering. 


2.0  HEAT  TRANSFER  ANALYSES 


2.1  THEORETICAL  ANALYSIS 


Heat  transfer  In  the  vehicle  and  from  the  vehicle  can  occur  in 

four  ways: 

(1)  Convection  from  a  surface  to  the  surrounding  fluids; 

(2)  Radiation  from  a  surface  directly  to  another  surface; 


(3)  Conduction  through  a  body  from  one  side  to  another; 


(4)  Convection  by  air  exchange  between  the  vehicle  and  the 
ambient  air  (ventilation). 

The  following  analysis  uses  the  first  three  heat  transfer 
mechanisms  to  predict  heat  flow  rates  and  the  effects  of  adding  insulation 
to  the  surfaces  of  the  vehicle.  As  ventilation  Is  essentially  Independent 
of  the  vehicle  Insulation,  Its  analysis  will  be  covered  In  a  separate 
section.  The  ventilatory  heat  loss  acts  In  parallel  with  the  other  three 
heat  loss  mechanisms  and  thus,  all  four  mechanisms  must  be  considered 
simultaneously  when  determlng  the  actual  effectiveness  of  an  Insulating 
cover. 

Dimensions  taken  from  a  test  vehicle  at  DREO  were  used  in  the 
Idealized  model  of  an  APC  shown  in  Figure  1. 

Convective  heat  transfer  coefficients  depend  upon  the  fluid 
properties.  As  these  properties  are  temperature  dependent.  It  Is  customary 
to  evaluate  them  at  a  temperature  which  Is  the  average  of  the  wall  and 
fluid  temperatures.  For  the  purpose  of  this  study,  the  outside  ambient  air 
temperature  Is  assumed  to  be  -40*C  while  the  Internal  air  temperature  Is 
assumed  to  be  15*C.  Fortunately,  the  physical  properties  of  air  (Appendix 
A)  are  only  mildly  dependent  upon  temperature  In  the  above  temperature 
range  and  even  sizable  deviations  from  the  assumed  temperatures  have  only 


small  effects  upon  the  magnitudes  of  the  heat  transfer  coefficients.  As 
the  internal  wall  temperatures  are  unknown,  it  is  initially  assumed  that 
all  the  walls  are  at  -12.5#C.  Using  this  temperature,  a  heat  transfer 
coefficient  may  be  calculated.  If  the  subsequent  evaluation  of  the 

predicted  wall  temperature  varies  by  more  than  5°C  from  this  value,  the 
heat  transfer  coefficients  will  be  re-evaluated  at  this  new  wall 
temperature  and  the  calculation  of  wall  temperatures  will  be  performed 
again:  This  iterative  procedure  may  be  repeated  as  often  as  necessary,  re¬ 
evaluating  the  heat  transfer  coefficients  using  the  latest  wall  temperature 
prediction. 

Several  dimensionless  groups  are  used  in  characterizing  heat 
transfer  and  fluid  flow  in  an  attempt  to  generalize  empirical  results.  The 
dimensionless  groups  of  interest  to  this  study  are: 

Reynolds  No.:  Re  *  UL^  Ratio  of  (Inertial  Forces)  to  (1) 

v  (Viscous  Forces). 


Nusselt  No.  :  Nu  =  hL  Ratio  of  (convective  heat  transfer  (2) 

IT"  rate  from  the  surface  to  the  fluid)  to 
(Conductive  heat  transfer  through 
the  liquid). 


Stanton  No.  :  St  =  h  Ratio  of  (Energy  transfer  rate  to  (3) 

Cpl)  the  fluid)  to  (Energy  flowrate  of  the 
H  fluid  due  to  the  fluid  velocity). 


Prandtl  No.  :  Pr  = 

Grashof  No.  :  Gr  = 


C_  v 

-E-  Ratio  of  (Molecular  diffuslvlty  of 

k  momentum)  to  (Molecular  dlffusivity 

of  heat). 


9pL3(T<-T  ) 


v2 


(4) 

(5) 


It  Is  convenient  to  analyze  the  problem  by  making  an  electrical 
analogy  where  the  temperature  difference  (which  is  the  driving  potential) 
Is  equivalent  to  the  voltage  drop  in  the  circuit;  the  heat  flow  rate  is 
equivalent  to  electrical  current;  and  the  thermal  heat  transfer 
coefficients  and  conductances  can  be  formulated  as  thermal  resistances 
corresponding  to  electrical  resistances  in  the  analogy. 

Figure  2  Is  a  graphical  representation  of  the  electrical  circuit 
which  is  equivalent  to  the  thermal  problem  of  a  heated  enclosure  which  has 
several  different  layers  of  thermal  resistance  between  the  interior  and  the 
exterior.  Each  large  "R"  represents  a  conductive  thermal  resistance;  each 
small  "r"  represents  a  radiative  thermal  resistance;  each  large  "R"' 
represents  an  equivalent  thermal  resistance  of  convective  and  radiative 
resistances  which  act  In  parallel. 

The  T^j  refer  to  surface  temperatures:  T^  are  the  surface 
temperatures  of  the  Interior  walls;  T^  are  the  surface  temperatures  of  the 


INTERIOR  AIR 
TEMPERATURE:  T0 


EXTERIOR  AIR 
TEMPERATURE:  T00 


Figure  2:  Electrical  Circuit  Analogy  of  the  Heat  Transfer  for 
the  Armoured  Personnel  Carrier. 
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The  electric  heaters  typically  produced  1300  [W]  each.  Each 
heater  had  an  electric  fan  which  produced  a  typical  mass  flowrate  of 
7.1xl0~2  [kg/s].  One  of  the  three  heaters  was  positioned  so  that  the  air 
was  blown  horizontally  while  the  other  two  were  positioned  so  that  the  air 
was  blown  upwards.  All  three  electric  heaters  in  each  vehicle  were  placed 
at  the.  front  of  the  crew  compartment  adjacent  to  the  vehicle  heater. 

The  thermistors  were  distributed  throughout  the  vehicle  so  as  to 
give  a  representative  temperature  distribution.  The  location  of  each 
thermistor  corresponded  as  closely  as  possible  to  the  same  location  in  each 
vehicle.  A  typical  thermistor  distribution  is  shown  in  Figure  3.  The  heat 
flow  discs  were  placed  on  the  floor  in  the  middle  of  the  cargo  space.  Both 
of  the  engine  compartment  doors  were  removed  to  allow  the  internal  air  to 
circulate  into  the  engine  compartment. 

Several  sets  of  measurements  were  taken,  each  measurement  period 
lasting  from  1800  hrs  to  0700  hrs  the  next  day.  In  case  1  of  the 
measurements,  neither  of  the  vehicles  had  additional  insulation.  This  was 
done  to  observe  any  differences  in  the  two  vehicles,  or  in  the  data 
collection  and  interpretation  technique.  In  cases  2a  and  2b,  the 
differences  between  the  heat  transfer  properties  of  an  uninsulated  APC  and 
an  APC  insulated  with  the  DRES  "Blanket"  are  shown.  In  case  3,  air  was 
allowed  to  Infiltrate  into  the  vehicle  through  the  hatches  to  observe 
ventilatory  heat  losses. 

The  data  were  analysed  by  assigning  a  volume  element  to  each  of 
the  thermistors  in  the  vehicle,  so  that  the  sum  of  all  the  volume  elements 
equalled  the  total  volume  of  the  APC.  Then,  using  these  volumes,  the 
volume  weighted  mean  temperature  of  the  interior  of  the  APC  was  calculated 
[7]  by: 


Ta  =  f  tiVv  (23) 

The  equivalent  thermal  resistance  (Req)  on  a  unit  area  basis  is  then 

calculated  by  dividing  the  difference  between  the  mean  temperature  and  the 
ambient  temperature  by  the  rate  of  heat  input: 


Req  =  A<Ta  “  TJ/Q  (24) 

where  the  total  vehicle  surface  area,  A,  is  assumed  to  be  42.125  m2.  Using 
this  parameter,  the  effectiveness  of  the  insulating  layer  can  be  assessed. 
Table  7  of  section  3.2  lists  tr.c  Req  of  the  vehicles  and  several  other 
Important  temperatures  found  in  the  experiments.  When  ventilation  of  the 
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2.2  ANALYSIS  OF  THE  EFFECT  OF  VENTILATION 


If  any  of  the  hatches  are  open,  or  other  openings  between  the 
interior  and  exterior  of  the  vehicle  exist,  thermal  energy  may  be  lost  by 
the  escape  of  warm  air  to  the  environment  which  is  then  replaced  by  cold 
ambient  air. 

If  it  is  assumed  that  the  rate  at  which  warm  air  escapes  from  the 
vehicle  is  q  [m3/s],  then  the  additional  vehicle  heating  required  to 
maintain  the  same  internal  temperature  is: 


Qv  -  P  Cpq(T0-Tj  (22) 


It  should  be  stressed  that  the  quantity  Qv  will  only  ensure  that 
the  mean  Internal  air  temperature  will  remain  unchanged.  Incomplete  mixing 
and  poor  circulation  will  undoubtebly  result  in  areas  of  cold  drafts  and 
hot  stagnant  pockets  of  air. 

The  only  current  practical  solution  to  this  heat  loss  problem  is 
to  minimize  the  air  exchange  with  the  environment.  The  lower  limit  of  this 
minimization  Is  the  rate  at  which  air  is  used  by  the  occupants  of  the 
vehicle  and  by  equipment  which  use  air,  such  as  combustion  heaters,  which 
are  operating  in  the  vehicle. 


2.3  EXPERIMENTAL  ANALYSIS 


Two  APC's  were  used  at  DREO  for  the  experimental  testing  of  the 
DRES  "Blanket".  Each  APC  was  made  "air-tight"  either  by  extensive  taping 
over  the  holes,  or  by  filling  the  holes  with  a  foam  insulating  compound. 
Missing  hatch-covers  were  replaced  with  a  6  mil  plastic  sheet  firmly 
secured  about  the  hatch  with  tape. 

After  the  APC's  were  made  as  air-tight  as  possible,  the  DRES 
Blanket  was  placed  on  one  of  the  vehicles.  Thermistors,  heat-flow  discs 
and  electrical  heaters  were  then  installed  inside  the  vehicles. 


f ABLE  4c:  Numerical  Values  of  Thermal  Resistances  for  Idealized 
APC:  Conductive  Heat  Transfer  for  Various  Insulating 
Configurations. 


I 


I 


Case  Description 


Case  1:  No  Insulating  cover. 

Case  2:  Idealized  ORES  "Blanket",  no  floor  insulation. 

Case  3:  Idealized  DRES  "Blanket"  plus  polyethylene  foam,  floor  insulation. 
Case  4:  Double  the  Insulation  thicknesses  of  Case  3. 

Case  5:  Quadruple  the  Insulation  thicknesses  of  Case  3. 
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TABLE  4a:  Initial  Numerical  Values  of  Thermal  Resistances  for  an 
Idealized  APC:  Internal  Convective  (R^),  Internal 
Conductive  (R^)  and  and  Exterior  Convective/Radiative 
(R14)  Heat  Transfer. 

{T0  =  15°C),  T„  =  -40 ®C,  Tm  *  -40 ®C,  V  *  20  m/s). 


Region 

Designator 

1 

R 

11 

(K/W) 

R 

12 

(K/W) 

R' 

14 

(K/W) 

Roof 

1 

3.1xl0“2 

4.5xl0-5 

2.2x10-** 

Side  Walls 

2,3 

5.0xl0-2 

8.8x10*5 

4.3x10-** 

End  Walls 

4,5 

8.1x10-2 

14.2x10*5 

9.5x10“** 

Floor 

6 

12.4x10-2 

4.5x10-5 

2.2x10'** 

TABLE  4b:  Numerical  Values  of  Thermal  Resistances  for  an  Idealized 
APC:  Internal,  Radiative  Heat  Transfer.  (Tm  *  15°C ) . 


\To 

From 

rij 

[k/W] 

j-1 

ja2 

ja3 

j=4 

ja5 

j*6 

Roof, 

1*1 

09 

0.1158 

0.1158 

0.2001 

0.2001 

0.0392 

Side  Wall  1, 

1*2 

0.1158 

09 

0.2305 

0.3660 

0.3660 

0.1158 

Side  Wall  2, 

1*3 

0.1158 

0.2305 

CD 

0.3660 

0.3660 

0.1158 

End  Wall  1, 

1*4 

0.2001 

0.3660 

0.3660 

09 

1.0813 

0.2001 

End  Wall  2, 

1*5 

0.2001 

0.3660 

0.3660 

1.0813 

09 

0.2001 

Floor, 

1*6 

0.0392 

0.1158 

0.1158 

0.2001 

0.2001 

09 

2.1.6  Numerical  Evaluation  of  Individual  Thermal  Resistances 


Tables  4a,  4b  and  4c  list  the  numerical  values  of  the  thermal 
resistances  of  the  model,  using  the  heat  transfer  coefficients  calculated 
in  Sections  2.1.1  to  2.1.3  with  the  appropriate  thermal  resistance  formulae 
and  the  radiant  thermal  resistances  of  Sections  2.1.4  and  2.1.5.  At  the 
exterior  surface,  the  convective  and  radiant  heat  transfer  act  in  parallel 
(assumption  2  of  Section  2.1.5).  Thus,  the  equivalent  thermal  resistance 
of  this  pair  nwy  be  calculated  as: 


R 


1 


1 

<W  +  h1Ai> 


(21) 


To  study  the  effects  of  adding  insulation  (varying  Ri3  of  Figure 
2),  five  different  Insulating  sample  cases  were  studied. 

Case  1:  No  additional  Insulation  (R^3  =  0  [K/W],  1*1,6); 

Case  2:  The  equivalent  of  the  DRES  "Blanket"  (R^3  =  0.424/A^  £, 

r63  *  0 


Case  3:  A  floor  insulation  equivalent  to  2  cm  of  closed  cell 
polyethylene  foam  (Rg3  =  0.441  5^-)  was  assumed  in 
addition  to  the  insulation  values  of  case  U 


Case  4:  Double  the  Insulating  values  of  Case  3; 
Case  5:  Quadruple  the  insulating  values  of  Case  3. 


14 


that  all  of  Its  Interior  walls  are  at  similar  temperatures,  and  these 
temperatures  must  not  be  greatly  different  from  the  surface  temperature  of 
the  bodies  Inside  the  vehicle.  (Note  that  the  surface  temperature  of  a 
clothed  soldier  nwy  be  significantly  different  from  skin  temperature. ) 


2.1.5  Radiant  Heat  Transfer  -  External 


In  general,  the  ground  and  sky  temperatures  can  be  significantly 
different  from  the  ambient  air  temperature  at  ground  level  [6].  The 
following  assumptions  were  made  to  simplify  the  analysis: 

1.  The  emlsslvlty  of  the  surroundings  Is  1.0; 

2.  The  temperature  of  the  surroundings  Is  the  same  as  the 
ambient  temperature. 

These  assumptions  may  be  refined  If  additional  information  Is  available. 
If  the  analysis  Is  to  be  applied  to  warm  weather  conditions  with  a 
significant  amount  of  solar  heating,  assumption  2  may  be  Invalid  and  a  more 
detailed  analysis  may  be  required. 

From  the  geometry  of  the  model,  the  vlewfactor  from  each  surface 
of  the  model  to  the  surroundings  Is  1.  Thus,  the  radiant  heat  transfer 
from  surface  1  to  the  surroundings  will  be  approximately: 


Q1  =  4A1cTm3(Ti  -  TJ  (19) 


and  the  thermal  resistance  for  the  surface  1  will  be: 


l 

r  *  ■  — ■  - 
1«  4A.j  crTm3 


(20) 


In  this  analysis,  the  radiant  heat  transfer  was  made  linear  In  temperature 
In  the  same  manner  as  In  Section  2.1.4.  As  a  first  order  approximation,  T_ 
Is  assumed  to  have  the  same  value  as  the  ambient  temperature  (-40°C). 


Calculated  vlewfactors  (Fjj)  are  given  In  Table  3.  When  these 
values  are  used  with  equation  (17 ) J  to  determine  the  radiant  thermal 
resistance.  It  Is  assumed  that  the  emlsslvlty  of  the  surfaces  are  all 
unity. 


TABLE  3:  Vlewfactors  for  Radiative  Heat  Transfer  Between 
Internal  Surfaces  of  an  Idealized  APC 


\to 

From 

Roof 

Side  Wall 

1 

Side  Wall 

2 

End  Wall 

1 

End  Wall 

2 

Floor 

Roof 

0 

0.1625 

0.1625 

0.0940 

0.0940 

0.4800 

Side  Wall 

1 

0.3185 

0 

0.1600 

0.1008 

0.1008 

0.3185 

Side  Wall 

2 

0.3185 

0.1600 

0 

0.1008 

0.1008 

0.3185 

End  Wall 

1 

0.2972 

0.1625 

0.1625 

0 

0.0550 

0.2972 

End  Wall 

2 

0.2972 

0.1625 

0.1625 

0.0550 

0 

0.2972 

Floor 

0.4800 

0.1625 

0.1625 

0.0940 

0.0940 

0 

It  should  be  noted  that  Internal  radiant  heat  transfer  can  be 
dominant  over  the  convective  heat  transfer  for  some  surfaces  If  the 
surfaces  have  significantly  differing  temperatures.  This  will  be 
especially  true  If  a  warm  body  Is  present  within  the  vehicle,  such  as  an 
engine  block  or  a  soldier.  If  the  surface  temperature  of  the  warm  body  Is 
approximately  the  same  temperature  as  the  ambient  air,  convective  heat 
transfer  will  be  negligible.  However,  If  the  vehicle  wall  temperatures  are 
significantly  different  from  the  air  temperature,  radiant  heat  transfer 
from  the  warm  body  will  occur  at  an  approximate  rate  of  1.5  W/m2  per  degree 
Celsius  of  temperature  difference  between  the  warm  body  and  the  walls. 
Thus,  for  an  effective  solution  to  the  problem  of  maintaining  a  warm  body 
Inside  the  vehicle  at  a  fixed  temperature,  the  vehicle  must  be  insulated  so 
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If  the  temperature  differences  between  the  surfaces  are  not  great,  then  the 
following  approximation  may  be  made: 


<T1H  •  v>  ■  4T»(T1  -  tj> 


(15) 


where 


Tn,  *  <T1  +  Tj  1/2.0 


(16) 


Thus  equation  (14)  becomes: 


0(  *  W  I  F,j(T(  -  Tj) 


(17) 


and  an  electrical  analogy  may  be  made  which  Is  similar  to  that  in 
convective/conductive  heat  transfer.  In  this  case,  Qj  is  composed  of 
several  heat  transfers  directly  to  the  other  walls,  so  the  radiative 
thermal  resistance  acts  In  parallel  to  the  conductive/convective  thermal 
resistance.  The  radiative  thermal  resistance  between  two  surfaces  1  and  j 
may  be  defined  by: 


rij '  <4VVF,j) 


(18) 


Unfortunately,  the  problem  Is  still  non-linear  due  to  the  Tm3  term.  To 
eliminate  the  non-linearity,  Tm  Is  Initially  assumed  to  be  equal  to  the 
Internal  air  temperature.  If  the  difference  between  any  wall  temperature 
and  the  Internal  air  temperature  Is  less  than  5*C  then  the  error  Incurred 
by  this  approximation  Is  less  than  5%.  For  applications  In  which  greater 
accuracy  Is  required,  Tm  may  then  be  re-evaluated  using  the  latest 
prediction  of  the  wall  temperatures  and  then  Iterating  the  solution  using 
the  Improved  values  of  the  radiative  thermal  resistance.  This  procedure 
will  also  be  used  for  external  radiative  heat  transfer  calculations. 


Table  2  lists  the  results  of  the  calculation  of  the  forced 
convective  heat  transfer  coefficient  for  the  outside  of  the  APC. 


TABLE  2:  Heat  Transfer  Coefficients  for  External,  Turbulent  Forced 
Convection,  Heat  Transfer  from  an  Idealized  APC. 

Free  Stream  Velocity,  V  =  20  m/s 


Region 

Area 

(m2) 

L 

(m) 

St 

h 

(W/m2K) 

Roof 

12.0 

4.2 

1.31x10*2 

3.76xl02 

Side  Wall 

6.1 

4.2 

1.31x10-2 

3.76x102 

End  Wall 

3.8 

2.9 

9.57xl0-3 

2.75x102 

Floor 

12.0 

4.25 

1.31x10-2 

3.76xl02 

2.1.4  Radiative  Heat  Transfer  -  Internal 


The  amount  of  radiant  energy  leaving  a  black  surface  (e*l)  Is 
given  by  [2,5]: 


Q,  -  A, off,*  -  J  F,jTj-> 


By  conservation  of  energy,  for  any  surface  1  the  sum  of  the  vlewfactors 
must  be  one.  Equation  13  therefore  reduces  to: 


q,  ■  vj 
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2.1.3  Convective  Heat  Transfer  -  External 


In  general,  it  cannot  be  expected  that  the  air-flow  over  the 
vehicle  will  be  small.  Indeed,  by  the  reasoning  put  forth  in  Section 
2.1.1,  if  the  air-speed  over  the  vehicle  Is  much  larger  than  1  [m/s]  (3.6 
[km/hr])  then  the  forced  convection  will  be  dominant.  Improper  use  of  a 
free  convection  analysis  severely  underestimates  the  heat  transfer  from  the 
outside  surfaces. 

In  the  analysis,  the  following  simplifying  assumption  Is  made: 
each  surface  is  considered  as  a  flat  plate  over  which  the  air  Is  blowing  at 
the  stated  free-stream  velocity.  This  assumption  will  over-estimate  the 
heat  transfer  coefficient  from  the  vehicle,  however,  as  the  outer  boundary 
layers  of  air  are  found  to  contribute  little  to  the  total  thermal 
resistance  of  the  vehicle,  the  errors  Induced  by  this  assumption  are 
thought  to  be  negligible. 

For  flow  along  a  flat  plate,  the  flow  begins  transition  at  a 
Reynolds  number  of  approximately  1x10$  and  becomes  fully  turbulent  when  the 
Reynolds  number  reaches  2.8xl06  [4].  If  a  wind  speed  of  20  [m/s]  (72 
[km/hr])  Is  present,  nearly  all  of  each  surface  (  88%)  will  be  subjected  to 
either  transitional  or  turbulent  flow.  For  the  purposes  of  this  study.  It 
will  be  assumed  that  each  surface  experiences  turbulent  flow  over  Its 
entire  area. 

Using  the  Reynolds  analogy,  the  local  Stanton  number  can  be 
related  to  the  Reynolds  number  by  [2]: 


1 


0.592(Rex)  5, 

5xl05<Rex<lxl07 

(11a) 

0.37(logRexr2'584, 

lxl07<Rex<lxl09 

(lib) 

The  mean  Stanton  number  can  be  determined  by  Integrating  the 
above  relations  over  the  length  of  the  surface.  For  convenience,  a  wind 
speed  of  20  [m/s]  was  chosen  which  produces  a  maximum  Reynolds  number  of 
8.9xlo6  which  Is  In  the  first  range  noted  above.  This  relationship  Is 
easily  Integrated  and  found  to  be: 


-0.2 


St  -  0.074  L(ReL) 


(12) 


TABLE  1:  Initial  Estimates  of  Heat  Transfer  Coefficients  for 
Internal,  Free  Convection,  Heat  Transfer  In  an 
Idealized  APC 


Region 

Area 

(m2) 

L 

(m) 

GrPr 

Flow 

Type 

Nu 

h 

(W/m2K) 

Roof 

9.8 

3.225 

3. 09x10 10 

Turbulent 

4.39xl02 

3.29 

Side  Wall 

5.0 

1.25 

8.81xl09 

Turbulent 

2.07xl02 

4.00 

End  Wall 

3.1 

1.25 

8.81xl09 

Turbulent 

2.07xl02 

1 

Floor 

9.8 

3.225 

3.09xl09 

Laminar 

1.13xl02 

0.85 
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vehicle.  If  It  Is  found  that  larger  air  velocities  over  the  Interior 
surfaces  are  present,  then  the  conclusions  drawn  from  the  following 
analysis  may  no  longer  be  valid. 

For  free  convection,  the  mean  Mussel t  number  has  been  found  to 
correlate  well  on  a  simple  function  of  the  Grashof  and  Prandtl  numbers  [2]: 


Nu  =  c(GrPr)n,  lO-^GrPr  <1013 


which,  when  evaluated,  gives  the  heat  transfer  coefficient: 


h  -  Nu  x  k 


The  constants  "c"  and  "n"  may  be  obtained  from  Table  B-l  for  free 
convection  from  vertical  surfaces  and  from  Table  B-2  for  free  convection 
from  horizontal  surfaces,  which  are  found  in  Appendix  B.  Table  1  gives  the 
results  of  the  analysis  for  the  internal  convective  heat  transfer 
coefficients. 


The  thermal  resistance  of  this  boundary  layer  of  air  can  be 
calculated  for  each  wall  from  the  equation: 

R  *  ~  (9) 


2.1.2  Conductive  Heat  Transfer  Through  the  Walls 


The  thermal  resistance  of  a  wall  (R)  having  a  known  thermal 
conductivity  (k),  thickness  { t )  and  area  (A)  is  given  by: 


r  * 


It  is  assumed  that  the  thermal  conductivity  of  the  aluminum  [3]  in  the 
walls,  floor  and  roof  of  the  APC  has  a  thermal  conductivity  of  204  [W/mK] 
and  that  all  of  these  surfaces  have  a  thickness  of  0.1  [»]. 


•  .•  •  ^ 
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A  computer  program  has  been  written  (Appendix  C)  which  solves 
this  set  of  simultaneous  equations  given  the  values  of  thermal  resistance 
and  ambient  temperatures  required  to  evaluate  the  elements  of  [A]  and  [B]. 
The  program  provides  the  total  heat  transfer  rate  for  the  enclosure  as  well 
as  the  Individual  heat  transfer  rates  through  each  surface.  Interior 
surface  temperatures  are  also  computed  and  presented  (Appendix  D). 

The  remainder  of  this  section  Is  devoted  to  the  evaluation  of  the 
heat  transfer  coefficients  and  the  thermal  resistances  for  the  model. 


2.1.1  Convective  Heat  Transfer  -  Internal 


Heat  transfer  by  convection  falls  within  two  categories:  Free 
Convection  and  Forced  Convection.  In  free  convection,  the  fluid  velocity 
generated  by  the  fluid  temperature  gradient  Is  dominant  over  any  Impressed 
fluid  velocity.  In  forced  convection,  the  opposite  Is  the  case. 

For  most  of  the  surfaces  In  the  vehicle,  any  impressed  velocity 
will  move  the  air  in  a  direction  which  Is  perpendicular  to  the  direction 
caused  by  natural  convection.  Thus,  the  resulting  heat  transfer  could 
conceivably  be  significantly  greater  than  that  predicted  by  either  analysis 
singly.  It  Is  difficult  to  determine  which  heat  transfer  mode  is  dominant 
or  even  significant  due  to  the  lack  of  data  on  air  velocities  inside  the 
vehicle.  It  has  been  observed  in  a  tentage  research  study  [1]  that  air 
velocities  across  the  floor  and  next  to  the  walls  are  typically  in  the 
range  of  0  to  0.5  [m/s]  for  Impressed  circulation  rates  of  approximately 
1  [m3/s].  A  measure  of  the  Importance  of  each  mode  of  convective  heat 
transfer  Is  given  by  the  ratio  of  the  Grashof  number  to  the  square  of  the 
Reynolds  number.  This  Is  Indicative  of  the  Importance  of  buoyancy  forces 

to  inertial  forces.  If  this  ratio  is  of  the  order  of  one,  then  both 

convection  modes  are  Important;  If  It  Is  significantly  greater  than  one, 
then  free  convection  Is  dominant.  The  magnitude  of  the  ratio  Is 

proportional  to  the  Inverse  of  the  square  of  the  free-stream  air  velocity. 
The  value  of  this  ratio  was  found  to  be  approximately  eight  when  an 
Impressed  air  velocity  of  0.5  [m/s]  was  assumed.  This  would  indicate  that 
both  forms,  natural  and  forced  convective  heat  transfer,  are  approximately 
of  equal  Importance  and  neither  one  should  be  disregarded.  Unfortunately, 
as  noted  above,  there  is  no  Information  on  Internal  APC  air  velocities  on 
which  to  make  valid  estimates  of  the  forced  convective  heat  transfer. 
Therefore,  for  the  purpose  of  this  study,  only  the  free  convection  will  be 
considered  with  the  understandl ng  that  it  may  underestimate  the  Internal 
heat  transfer.  When  the  required  data  on  Internal  air  velocities  are 
available.  It  Is  recommended  that  this  portion  of  the  analysis  be  examined 
to  determine  the  error  Introduced  by  considering  only  the  free  convection 
and  neglecting  the  forced  convection  on  the  Internal  surfaces  of  the 


outside  walls  and  the  inner-most  surfaces  of  any  insulating  cover;  are 
the  outside  surface  temperatures  of  any  Insulating  cover. 

Using  conventional  circuit  analysis  techniques ,  the  problem  Is 
reduced  to  solving  a  set  of  simultaneous  equations.  The  Inside  wall 
temperatures  of  matrix  [T]  are  the  primary  unknowns  which  can  be  expressed 
in  terms  of  the  known  thermal  resistances  of  matrix  [A]  and  the  known 
Internal  and  external  temperatures  of  matrix  [B]. 

At  each  of  the  6  nodes  (11,21,31,41,51,61),  which  correspond  to 
six  interior  walls,  a  heat  flow  balance  to  other  nodes  Is  performed: 


6 

%1  =  ^1-  +  ji^lj 


The  heat  flows  are  then  replaced  by  the  appropriate  temperature  differences 
and  thermal  resistances  between  nodes  according  to  the  relationship: 


q 


,4L 

R 


Upon  separation  of  the  Individual  wall  temperatures  and  simplifying,  a  set 
of  simultaneous  equations  Is  found: 


[A][T]  =  [B] 


(6a) 


where. 


A  ( 1 ,  j )  =  + 


'll  {R12  +  Ri3  +R  * 14J 


+  ?  1  ■  j 

*  r1k 


=  A 

r 


B  ( 1 )  =  -  (- 


1j 


'll  (R12  +  R13  +  R14} 


-  j 

(6d) 

*  3 

(6c) 

(6d) 

Having  obtained  the  temperatures  of  [T],  It  Is  a  simple  procedure  to  obtain 
the  various  heat  flow  components.  The  procedure  Is  completely  analogous  to 
determining  the  unknown  currents  In  an  electrical  circuit  analysis. 


TOP  VIEW 


Typical  Experimental  Layout  of  Sensors  and  Equipment 
in  the  Armoured  Personnel  Carriers. 
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vehicles  occurs,  Req  is  a  function  of  both  the  ventilation  rate  as  well  as 
the  rate  at  which  heat  is  lost  by  convection,  conduction  and  radiation. 


3.0  RESULTS 


Table  5  contains  the  analytical  estimates  of  the  thermal  energy 
requirements  to  maintain  a  15°C  internal  temperature  with  an  external 
ambient  temperature  of  -40#C.  It  is  assumed  that  no  ventilation  occurs. 
These  values  are  taken  or  derived  from  the  data  listed  in  Appendix  D.  The 
"Test  Number"  indicated  with  the  data  in  Appendix  D  is  chosen  such  that  the 
first  digit  corresponds  to  the  case  number  and  the  second  digit  indicates 
the  iteration  number. 

It  was  found  that  the  value  of  the  wall  temperature  initially 
assumed  was  often  significantly  different  from  the  predicted  wall 
temperature.  Thus,  the  heat  transfer  coefficients  were  recalculated  using 
improved  wall  temperature  estimate  and  the  analysis  was  repeated  until  the 
difference  between  two  successive  predictions  of  wall  temperatures  were 
less  than  1"C.  This  difference  results  in  less  than  a  2%  error  in  the 
calculation  of  the  heat  transfer  coefficients.  Three  iterations  were 
usually  required  to  achieve  this  precision. 

TABLE  5:  Predicted  Values  of  Heating  Requirement(Q),  Mean  Wall 
Temperature  (T_ )  and  Equivalent  Thermal  Resistance  (R) 
of  an  Idealized  APC  for  each  Insulating  Case  defined 
in  Table  4c. 

(T0  =  15  *C,  T.  =  —40 °C ) . 


Case 

Q 

CM] 

Ts 

ra 

R 

£5] 

w 

1 

9515 

-39.3 

0.21 

2 

4174 

-20.7 

0.47 

3 

2866 

-  9.7 

0.69 

4 

1777 

-  2.4 

1.11 

5 

1030 

3.6 

1.91 

•  ■  '.X-''  .S"V-V  \ 
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For  each  Insulating  case  of  Table  4c,  a  mean  thermal  resistance  Is 
calculated  for  the  entire  vehicle  by  dividing  the  difference  between  the 
internal  and  ambient  air  temperatures  by  the  computed  heat  loss.  Also, an 
area-weighted  mean  Interior  surface  temperature  (T$)  Is  calculated. 

Examination  of  the  results  indicates  that  merely  considering  the 
thermal  energy  requirements  to  maintain  the  difference  between  the  Internal 
and  ambient  air  temperatures  will  given  an  Incomplete  view  of  the  heat 
transfer  problem.  The  importance  of  the  interior  wall  temperature  can  be 
studied  with  reference  to  the  five  cases  defined  earlier. 

The  analysis  of  Case  1  (uninsulated  vehicle)  Indicates  that  It 
requires  9515  [W]  to  maintain  the  internal  air  of  a  temperature  55*C  higher 
than  the  ambient  air.  However,  most  of  the  vehicle's  thermal  Insulation 
results  from  the  Internal  air  boundary  layers  which  form  on  the  vehicle 
walls.  This  results  In  a  mean  Internal  wall  surface  temperature  of 
-39.3*C  which  Is  approximately  equal  to  the  ambient  temperature  (-40*0. 
If  a  body  which  has  a  surface  temperature  which  Is  conparable  to  the 
Internal  air  temperature  Is  placed  Into  the  vehicle,  the  convective  heat 
transfer  from  that  body  will  be  negligible.  However,  thermal  energy  would 
be  lost  from  the  body  by  radiative  heat  transfer  to  the  walls  at  an  Initial 
rate  of  222  [W/m2].  As  a  guide,  the  heat  production  of  a  typical.  Inactive 
person  Is  approximately  45  [W/m2].  Thus,  It  Is  apparent  that  the  Internal 
wall  surface  temperature  should  be  approximately  equal  to  the  warm  body 
surface  temperature  to  minimize  radiative  heat  ransfer  from  the  body. 

In  the  second  case,  (modeling  the  addition  of  the  DRES 
"Blanket"),  the  heating  requirements  are  reduced  44%  from  the  uninsulated 
case.  The  highest  heat  loss  area  is  now  the  floor,  which  has  no  additional 
Insulation.  The  temperature  of  the  floor  (-39.4*0  Is  still  very  close  to 
the  ambient  air  temperature  compared  with  the  other  Interior  surface 
temperatures  (-14*C).  The  Initial  radiative  heat  loss  from  the  body  In  the 
vehicle  would  be  reduced  to  approximately  161  [W/m2].  By  adding  an 
Insulating  covering  to  the  floor  (the  equivalent  of  2.2  cm  of  a  closed 
cell,  polyethylene  foam  pad)  the  heat  loss  would  be  reduced  a  further  14% 
from  the  uninsulated  case,  but  the  floor  temperature  would  be  Increased 
significantly  to  -12°C  and  the  other  wall  temperatures  would  also  Increase 
substantially  to  -10*C.  This  step  would  reduce  the  initial  radiative  heat 
loss  from  the  body  to  118  [W/m2]. 

Insulating  cases  4  and  5  demonstrate  the  effects  of  doubling  and 
quadrupling  the  thickness  of  Insulation  of  case  3.  Case  5,  which  Is  the 
equivalent  of  6  cm  of  the  rayon  Insulation  of  the  DRES  "Blanket"  and  8  cm 
of  the  closed  cell  foam  pads,  allows  the  vehicle  to  maintain  the  55*C 
temperature  difference  with  only  11%  of  the  energy  required  for  uninsulated 
vehicle.  The  mean  wall  temperature  In  the  vehicle  for  this  case  Is 
approximately  3.6*C,  or  43.6*C  above  ambient.  Thus,  the  steady-state 
temperature  of  a  non-heat-producing  body  placed  In  the  vehicle  would  be 
somewhere  between  3.6*C  and  15*C.  For  people,  this  higher  wall  temperature 
would  mean  less  heat  loss  by  radiation  and  hence  a  more  comfortable 


23 


environment.  The  Initial  radiative  heat  transfer  from  the  boc|y  would  be  87 
and  58  [W/m2]  for  cases  4  and  5  respectively. 

Table  6  gives  examples  of  the  additional  heating  requirements  If 
there  were  to  be  an  air  exchange  between  the  Interior  of  the  vehicle  and 
the  surroundings.  The  ventilation  rate  will  be  a  complex  function  of  the 
number  and  size  of  openings  In  the  vehicle,  the  wind  velocity,  and  the 
temperature  difference. 


TA8LE  6:  Predicted  Heat  Loss  Rates  for  Hypothetical 
Air  Infiltration  Rates. 

(Tq  *  15 *C,  T.  =  -40 *C ) 


Infiltration  Rate 

Rate  of  Heat  Loss 

0n3/s] 

[W] 

0.01 

837.9 

0.10 

8379.2 

1.00 

83791.6 

Since  no  data  Is  currently  available  on  the  ventilation  rates,  a 
typical  value  was  assumed  In  order  to  demonstrate  the  significance  of  this 

form  of  heat  loss.  A  minimum  ventilation  rate  was  derived  by  assuming  that 

each  person  In  the  vehicle  requires  1  1/s  of  fresh  air  to  keep  the  carbon 
dioxide  levels  In  the  vehicle  below  1%,  and  by  assuming  that  there  were  ten 
people  In  the  vehicle.  Thus,  the  minimum  ventilation  rate  would  typically 
be  10  1/s  or  10“2  [m3/s].  From  Table  6,  It  is  seen  that  this  requires  an 
additional  838  W  of  thermal  energy  be  supplied  to  maintain  the  same  mean 
air  temperature  as  when  no  ventilation  occurred.  If  the  hatches  were  open. 
It  Is  expected  that  the  ventilation  rate  could  be  much  higher  than  the 
value  used  In  this  example.  Again,  It  must  be  emphasized  that  this 
additional  energy  Input  would  only  preserve  the  mean  air  temperature  and  It 
can  be  expected  that  regions  of  very  cold  and  very  warm  air  would  be 

present  In  the  vehicle. 

The  results  of  the  experimental  investigation  of  the 

effectiveness  of  the  ORES  "Blanket"  are  given  in  Table  7.  In  most  cases, 
the  weather  was  clear  with  no  wind.  As  previously  noted,  several  different 
"insulating"  conditions  were  Investigated  and  each  variation  Is  noted  In 
the  table. 


The  first  set  of  measurements  Indicate  that  there  are  only  minor 
differences  In  the  heat  transfer  properties  of  each  of  the  uninsulated 
vehicles.  The  calculated  value  of  the  equivalent  thermal  resistance  (Reg) 
Indicate  a  3.6*  difference  between  vehicles  and  It  Is  assumed  that  the  data 
analysis  technique  Is  appropriate. 
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The  value  of  Req  for  the  uninsulated  vehicle  (0.235[m2K/W])  Is  In  close 
agreement  with  that  predicted  by  the  analytical  model  (0. 2lCm2K/W3 ) - 


The  addition  of  the  DRES  "Blanket"  (cases  2a  and  2b)  Is  seen  to 
result  In  an  Increase  of  the  Rea  of  an  uninsulated  APC  by  approximately  24% 
wlth  .no  wind.  The  values  Of  Rea  (0.2908  and  0.3119  [m2K/W]  are 
significantly  less  than  the  predicted  value  (0.47  [m^/W]). 


Air  Infiltration  between  the  vehicle  and  the  "Blanket"  may  be 
responsible  for  this  discrepancy,  however,  air  movement  In  the  vehicle  and 
ventilation  may  also  be  responsible.  As  previously  noted,  air  movement  In 
the  vehicle  would  cause  heat  to  be  transferred  to  the  walls  at  a  greater 
rate  than  that  predicted  by  a  natural  convective  heat  transfer  analaysls. 


Comparison  of  the  data  for  heat  loss  through  the  floor  (Table  7) 
with  the  predicted  heat  loss  (Appendix  D)  on  a  per  unit  area  per  temper¬ 
ature  difference  gives  varying  agreement.  For  the  uninsulated  vehicles, 
the  average  heat  loss  through  the  floor  was  measured  to  be  0.96  [W/m3<] 
while  the  predicted  heat  loss  was  1.66  [W/m^].  For  the  vehicle  insulated 
with  the  DRES  "Blanket",  the  measured  heat  loss  through  the  floor  was 
measured  to  be  1.47  [W/m^]  while  the  predicted  heat  loss  was  4.17  [W/m?K]. 
The  discrepancy  between  the  measured  floor  heat  loss  and  the  predicted 
floor  heat  loss  is  probably  coupled  to  the  discrepancy  between  the 
theoretical  evaluation  of  the  heat  transfer  coefficient  and  the  Req  for  the 

vehicle  with  the  DRES  "Blanket".  From  test  case  2  (Appendix  D)  73*  of  the 
heat  transferred  to  the  floor  Is  by  radiation  from  the  other  walls.  If  the 
other  wall  temperatures  are  lower  than  expected  due  to  an  unexpectedly  low 
value  of  Req ,  then  the  radiant  heat  transfer  to  the  floor  may  be 

significantly  less  than  predicted. 


In  Case  3,  the  effect  of  ventilation  is  apparent.  Neither  APC's 
have  any  additional  Insulation.  The  first  APC  has  both  the  Crew 
Commander's  and  the  Driver's  hatches  open,  while  the  second  APC  has  only 
the  Driver's  hatch  open.  The  results  In  Table  7  Indicate  that  opening  the 
Driver's  hatch  causes  an  apparent  decrease  In  the  Req  of  an  uninsulated 

vehicle  of  17%.  Opening  both  the  Crew  Commander's  and  the  Driver's  hatches 
Is  seen  to  reduce  the  apparent  Req  of  an  uninsulated  vehicle  by  37%. 

In  Case  3,  the  wind  speed  was  negligible  and  thus  the  ventilation 
rate  resulted  from  the  bouyancy  of  the  Internal  air  only.  In  general.  If  a 
wind  were  present.  It  Is  expected  that  the  ventilation  rate  could  be 
substantially  Increased,  greatly  Increasing  the  heat  transfer  from  the 
vehicle.  Also,  as  the  Insulation  of  the  vehicle  is  Increased,  and  the  heat 
loss  through  the  walls  Is  reduced,  the  importance  of  ventilation  rate  and 
Its  associated  heat  loss  become  proportionally  larger. 


Results  Derived  from  Experimental  Measurements  of  Temperature  In  Two  APC 
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These  percentage  differences  are  calculated  based  on  the  equivalent  thermal  resistance  of  the 
correspondl ng  APC  relative  to  the  mean  equivalent  thermal  resistance  of  the  uninsulated  APC's 
measured  In  experiments  1,  2a,  2b  (Req  =  0.2352  m2K/W). 


4.0  CONCLUSION 


The  final  choice  of  the  amount  of  insulation  to  be  used  and  how 
it  is  to  be  distributed  on  the  vehicle  will  depend  upon  space  limitations, 
available  thermal  energy  and  the  desired  steac(y-state  temperature  of  bodies 
inside  the  vehicle.  Due  to  the  transfer  of  thermal  energy  by  both 
conduction  and  radiation,  it  is  apparent  that  the  insulation  should  be 
distributed  so  as  to  give  a  uniform  internal  wall  temperature  and  that  this 
interior  wall  temperature  should  be  as  close  to  the  internal  air 
temperature  as  possible. 

If  ventilation  occurs,  then  the  value  of  the  insulation  may 
quickly  become  of  secondary  Importance,  depending  upon  the  ventilation 
rate.  Even  at  a  minimum  ventilation  rate  of  10  1/s,  and  with  the 
insulation  of  case  5,  the  heating  requirements  for  the  ventilation  and 
conductive/radiative  heat  losses  are  of  the  same  order  of  magnitude, 
approximately  1000  W  each. 

Experimental  measurements  of  the  internal  temperatures  of  two 
APC's,  one  without  insulation  and  one  with  the  ORES  "Blanket",  Indicate 
that  this  covering  Increases  the  thermal  resistance  of  an  uninsulated 
vehicle  from  0.2352  to  0.3013  [m^/W].  This  would  allow  a  22%  reduction  in 
energy  Input  from  9851  [WJ  to  7690  [W]  while  still  maintaining  the  55°C 
temperature  difference.  Unfortunately,  this  may  be  of  little  benefit  if 
the  radiative  heat  transfer  from  a  boc(y  in  the  vehicle  Is  large  due  to  the 
low  floor  temperature.  From  the  analytical  stu<ty,  the  performance  of  the 
DRES  "Blanket''  can  be  significantly  enhanced  if  the  floor  Is  insulated  as 
well . 

The  discrepancies  between  Req  predicted  by  the  model  and  the  Re_ 
found  by  experimentation  for  the  Insulated  vehicle  are  significant.  When 
data  on  the  airflow  pattern  in  a  vehicle  are  available,  the  validity  of  the 
internal  heat  transfer  coefficients  should  be  determined  and  corrected  if 
necessary.  Some  of  the  discrepancy  may  be  In  the  manner  in  which  the 
insulating  covering  is  put  on  the  vehicle.  A  poor  fit  may  allow  air  to 
circulate  between  the  vehicle  and  the  covering,  reducing  the  effectiveness 
of  the  covering  as  a  thermal  Insulator.  This  effect  would  be  even  more 
significant  if  the  vehicle  were  to  be  exposed  to  a  wind.  Portions  of  the 
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vehicle  not  covered  by  the  blanket,  the  hatch  covers  for  example,  may  be 
sources  of  high  heat  loss. 

It  should  be  pointed  out  that  If  the  discrepancy  arises  from  the 
evaluation  of  the  Internal  heat  transfer  coefficients,  the  Importance  of 
the  error  decreases  as  more  Insulation  Is  added  to  the  vehicle.  If  the 
error  Is  a  result  of  air  Infiltration,  the  error  will  increase  with 
Increasing  insulation  and  It  may  vary  with  ambient  conditions.  If  the 
discrepancy  Is  a  result  of  a  poorly  fitted  covering  or  exposed  surfaces, 
little  change  in  the  thermal  performance  of  a  covering  should  be  observed 
as  the  insulating  value  of  the  covering  is  increased. 
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APPENDIX  A 

PHYSICAL  PROPERTIES  OF  AIR 


TABLE  Al:  Physical  Properties  of  Air  at  Standard 

Atmospheric  Pressure.  Derived  from  [8,9]. 


perature 

T 

C°C] 

Density 

p 

[kg/m3] 

Kinematic 

Viscosity 

[m2/s] 

Specific 

Heat 

[J/i?g  K] 

Thermal 

Conductivity 

K 

[W/mK] 

Prandtl 

Number 

Pr 

-40 

1.515 

0.98  x  10~5 

1005.6 

0.02086 

0.728 

-20 

1.395 

1.15 

1005.3 

0.02252 

0.721 

0 

1.293 

1.32 

1005.5 

0.02411 

0.716 

10 

1.248 

1.41 

1005.6 

0.02490 

0.713 

20 

1.205 

1.50 

1005.7 

0.02570 

0.710 

30 

1.165 

1.60 

1005.9 

0.02648 

0.707 

40 

1.128 

1.68 

1006.6 

0.02724 

0.705 
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APPENDIX  B 


EMPIRICAL  CONSTANTS  FOR  CONVECTIVE  HEAT  TRANSFER 


TABLE  Bl:  Constants  c  and  n  of  Eq.  (7)  for  Free  Convection 
from  a  Vertical  Plate  or  Cylinder  at  Uniform 
Temperature  [2]. 


Type  of  Flow 

Range  of  Gr|_Pr 

c 

n 

Laminar 

10**  to  109 

0.59 

1/4 

Turbulent 

109  to  1010 

0.10 

1/3 

TABLE  B2:  Constants  c  and  n  of  Eq.  (7)  for  Free  Convection  from  a 
Horizontal  Plate  at  Uniform  Temperature  [2]. 


Type  of  Flow 

Orientation  of  Plate* 

Range  of  GrLPr 

c 

n 

Laminar 

Upper  surface  heated 
or 

Lower  surface  cooled 

105  t0  2xl07 

0.54 

1/4 

Turbulent 

Upper  surface  heated 
or 

Lower  surface  cooled 

2x10 7  to  3xl010 

0.14 

1/3 

Laminar 

Lower  surface  heated 
or 

Upper  surface  cooled 

3><10 5  to  3*1010 

0.27 

1/4 

*  Action  of  the  fluid  on  the  plate. 


33 


APPENDIX  C 


PROGRAM  LISTINGS 


34 


1  1  C 

1  2  C 

1  3  PROGRAM  HTCALC 

1  4  C 

1  5  C 

1  6  C  THIS  PROGRAM  IS  DESIGNED  TO  SOLVE  THE  SYSTEM  OF  LINEAR, 

1  7  C  SIMULTANEOUS  HEAT  TRANSFER  EQUATIONS  BY  USING  A  LU  DECOMPOSITION, 

1  8  C  BACK-SUBSTITUTION  AND  AN  IMPROVEMENT  ROUTINE.  IT  IS  CURRENTLY 

1  9  C  SET  UP  TO  HANDLE  A  SET  OF  AS  MANY  AS  50  LINEAR  EQUATIONS 

1  10  C  WITH  50  UNKNOWN  TEMPERATURES.  DERIVATION  OF  THE  ELEMENTS 

1  11  C  OF  THE  MATRICIES  'A', 'TIN'  AND  'TOUT'  ARE  AS  GIVEN  IN  THE  TEXT. 


1 

12 

C 

THIS  PROGRAM  MAKES  USE  OF  THE  FOLLOWING  SUBPROGAMS: 

1 

13 

C 

DCOMS  - 

A  SUBROUTINE  FOR  THE  MATRIX  REDUCTION; 

1 

14 

C 

SOLVE  - 

A  SUBROUTINE  FOR  THE  BACK -SUBSTITUTION; 

1 

15 

C 

IMPRV  - 

A  SUBROUTINE  FOR  AN  ITERATIVE  IMPROVEMENT  SCHEME; 

1 

16 

C 

DTPROD- 

A  FUNCTION  WHICH  DETERMINES  THE  DOT  PRODUCT  OF 

1 

17 

C 

TWO 

VECTORS. 

1 

18 

C 

1 

19 

C 

1 

20 

C 

PARAMETERS: 

1 

21 

C 

1 

22 

C 

1 

23 

C 

A 

- 

A  2-D  ARRAY  WHICH  CONTAINS  THE  COEFFICIENTS 

1 

24 

C 

CORRESPONDING  TO  THE  UNKNOWNS  IN  THE  VECTOR  T. 

1 

25 

C 

THESE  COEFFICIENTS  ARE  COMPOSED  OF  THE  THERMAL 

1 

26 

C 

RESISTANCES  OF  THE  SYSTEM. 

1 

27 

C 

1 

28 

C 

AREA 

A  CHARACTER  ARRAY  WHOSE  ELEMENTS  ARE  SURFACE 

1 

29 

C 

NAMES  WHICH  CORRESPOND  TO  THE  SURFACE  DESIGNATING 

1 

30 

C 

ELEMENT  NUMBER. 

1 

31 

C 

1 

32 

C 

B 

- 

A  REAL  VECTOR  CONTAINING  THE  RIGHT  HAND  SIDE 

1 

33 

C 

OF  THE  SYSTEM  OF  EQUATIONS. 

1 

34 

C 

1 

35 

C 

IND 

- 

AN  INTEGER  INDICATING  WHETHER  OR  NOT  "A"  IS 

1 

36 

C 

SINGULAR.  IF  IND  =  -1  "A"  IS  SINGULAR  AND 

1 

37 

C 

YOU'VE  GOT  PROBLEMS.  IF  IND  =  0  THEN  "A" 

1 

38 

c 

IS  NONSINGULAR  AND  EVERYTHING  SHOULD  PROCEED. 

1 

39 

c 

1 

40 

c 

N 

- 

AN  INTEGER  INDICATING  THE  NUMBER  OF  UNKNOWN 

1 

41 

c 

SURFACE  TEMPERATURES  FOR  WHICH  THE  PROGRAM 

1 

42 

c 

WILL  SOLVE.  IT  IS  NECESSARY  THAT  N>*NDIM 

1 

43 

c 

SINCE  IF  N<NDIM  THE  PROBLEM  IS  INDETERMINANT. 

1 

44 

c 

1 

45 

c 

NDIN 

1  - 

AN  INTEGER  INDICATING  THE  NUMBER  OF  LINEAR 

1 

46 

c 

EQUATIONS  IN  THE  MATRIX  'A'. 

1 

47 

c 

1 

48 

c 

RCON 

1  - 

A  MATRIX  CONTAINING  THE  VALUES  OF  THE  CONDUCTIVE 

1 

49 

c 

AND  CONVECTIVE  THERMAL  RESISTANCES  AS  DEFINED 

1 

50 

c 

IN  THE  TEXT. 

1 

51 

c 

1 

52 

c 

RRAC 

1  - 

A  MATRIX  CONTAINING  THE  VALUES  OF  THE  RADIATIVE 

1 

53 

c 

THERMAL  RESISTANCES  AS  DEFINED  IN  THE  TEXT. 

1 

54 

c 

1  55  C  T  -A  REAL  VECTOR  WHICH  WILL  ULTIMATELY  CONTAIN 

1  56  C  THE  DESIRED  TEMPERATURES  (A*T=B) 

1  57  C 

I  58  C  TESTNO-  AN  INTEGER  IDENTIFYING  THE  TEST  NUMBER 

1  59  C 

1  60  C  TIN  -  THE  MEAN,  INTERNAL  AIR  TEMPERATURE  OF  THE  VEHICLE 

1  61  C 

1  62  C  TOUT  -  THE  EXTERIOR  AMBIENT  AIR  TEMPERATURE 

1  63  C 

1  64  C 

1  65  C 

1  66  C  DATA  INPUT  SCHEDULE: 

1  67  C 

1  68  C  RECORD  1;  TESTNO, NDIM,N 

1  69  C  RECORD  2-7;  AREA(I), 1=1,6 

1  70  C  RECORD  8-11;  (RCON( I , J) , 1=1 ,6) , J=l,4 

1  71  C  RECORD  12-16;  (RRAD(I,J)  ,1=1,5), >1+1,6 

1  72  C  RECORD  17;  TIN, TOUT 

1  73  C 

1  74  C 

1  75  REAL  A(50,50) ,B(50) ,T(50) ,D(50) ,Z(50) 

1  76  REAL  LU(50,50) ,R(50) ,RC0N(50,50) ,RRAD(50,50) , TIN, TOUT 

1  77  REAL  QC0N(50) ,QRAD(50,50) ,QT0T,RESIST(50) 

1  78  CHARACTER*12 ,  AREA(6) 

1  79  INTEGER  TESTNO, I, J,K,NDIM,N 

1  80  DIMENSION  NPIV(50) 

1  81  C 

1  82  READ(105,*)  TESTNO, NDIM,N 

1  83  C 

1  84  WRITE(*,5)  TESTNO 

1  85  5  FORMAT ( '  ', 'HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED' 

1  86  +'  APC  MODEL',/,'  TEST  NUMBER:  ',110,//) 

1  87  C 

1  88  DO  10  1=1,6 

1  89  READ(105,*)  AREA(I) 

1  90  10  CONTINUE 

1  91  C 

1  92  READ(105,*)  ( (RCON( I , J) , 1*1,6) , >1,4) 

1  93  C 

1  94  WRITE ( * , 15 )  ( ( AREA( I) ,RCON( I ,3) ), 1=1 ,6) 

1  95  15  FORMAT ( '  ADDED  INSULATION  VALUES  (K/W) ’ ,/, 

1  96  +6( '  R( ' ,A12, ' )"  ' , F10 . 6 , / ) ) 

1  97  C 

1  98  DO  30  1=1,5 

1  99  READ(105,*)  (RRAD( I ,K ) ,K=I+1 , 6) 

1  100  C 

1  101  DO  20  >1,6 

1  102  C 

1  103  IF  (I  .EQ.  J)  THEN 

1  104  C 

1  105  C  SET  RRAD(L,L)=1E99  WHICH  IS  INFINITY.  IT  IS  ACTUALLY 

1  106  C  INFINITE  FOR  NON-CONCAVE  SURFACES  AS  THE  VIEWFACTOR  FOR 

1  107  C  NON-CONCAVE  SURFACES  IS  ZERO. 

1  108  C 


v'  i  ■  *  ,  » ,  »  i 


ii  0  .  |.ii  v  i  j-n  tii  hip  ij  j 


'.v.’7.:v  v  v?  v 
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1 

109 

RRAD( I, J)  =  1.E99 

1 

no 

C 

1 

in 

ELSE 

1 

112 

RRAD( J, I)  =  RRAD(I.J) 

1 

113 

END  IF 

1 

114 

C 

1 

115 

20 

CONTINUE 

1 

116 

-  30 

CONTINUE 

1 

117 

C 

1 

118 

RRAD{6,6)  =  1.E99 

1 

119 

C 

1 

120 

READ(105 ,*)  TIN, TOUT 

1 

121 

C 

1 

122 

C 

NOW  COMPOSE  THE  ELEMENTS  OF  THE  MATRICES  'A1  AND  1 B ' 

1 

123 

C 

1 

124 

DO  70  1=1,6 

1 

125 

DO  60  0=1,6 

1 

126 

C 

1 

127 

IF  (I  .EQ.  J)  THEN 

1 

128 

C 

1 

129 

RES  1ST ( I)  =  0.0 

1 

130 

DO  40  K=2,4 

1 

131 

RESIST ( I)  =  RESIST ( I )  +  RCON(I,K) 

1 

132 

40 

CONTINUE 

1 

133 

A(I,J)  =  -(l./RCQN( 1,1)  +  1 ./RESIST ( I ) ) 

1 

134 

C 

1 

135 

DO  50  K=l,6 

1 

136 

A( I, J)=A( I, J)-1./RRAD( I,K) 

1 

137 

50 

CONTINUE 

1 

138 

C 

1 

139 

ELSE 

1 

140 

A(  I , J)  =  1 ./RRAD( I , J) 

1 

141 

A(J,I)  =  A(I,J) 

1 

142 

END  IF 

1 

143 

C 

1 

144 

60 

CONTINUE 

1 

145 

70 

CONTINUE 

1 

146 

C 

1 

147 

DO  90  1=1,6 

1 

148 

B( I )  =  -TIN/RCON( 1,1) -TOUT/RESIST ( I ) 

1 

149 

90 

CONTINUE 

1 

150 

C 

1 

151 

DO  110  1=1 ,NDIM 

1 

152 

DO  100  J=1,N 

1 

153 

LU( I , J)=A( I , J) 

1 

154 

100 

CONTINUE 

1 

155 

110 

CONTINUE 

1 

156 

C 

1 

157 

C 

NOW  PERFORM  THE  GAUSSIAN  ELIMINATION  (LU  DECOMPOSITION) 

1 

158 

C 

1 

159 

CALL  DCOMS(NDIM,N,LU,NPIV,D, IND) 

1 

160 

C 

1 

161 

C 

CHECK  TO  SEE  IF  THE  MATRIX  IS  SINGULAR 

1 

162 

C 

50 


HEA  r  TRANoFb-R  C ALcULA  f  i UN3  FOR  ttN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  22 


ADDED  INSULATION  VALUES  (K/W) 

K  ( RiJOF  —  ,  035290 

R(SIDE  WAlL  1  )=  .069430 

R03IDE  WALL  2  )=  .069430 

R  >  END  WALL  1  )=  .111450 

R (END  WALL  2  )=  .111450 

R (FLOOR  ) =  .000000 

CONDITION  NUMBER  =  1.693 

INTERNAL  WALL  SURFACE  TEMPERATURES  <C) 

T ( ROOF  )=  -12.39 

T  (SIDE  WALL  1  5=  -14.30 

T  (SIDE  WALL  2  )=  -14.30 

T (END  WALL  1  )=  -14.09 

F (END  WALL  2  >=  -14.09 

T  (FLOOR  )=  -39.41 

HEAT  LOSS  THROUGH  THE  WALLS  <W> 

Q ( ROOF  ) =  762.53 

0 < S IDE  WALL  1  >=  367.46 

GKSIDE  WALL  2  >=  367.46 

0 (END  WALL  1  )=  230.19 

O': END  WALL  2  )=  230.19 

O'.  FLOOR  )=  2216.35 

RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  <W> 


OR AD ( ROOF 

, SIDE  WALL  1 

>  = 

12. 

16 

OR AD ( ROOF 

.SIDE  WALL  2 

)  = 

12. 

16 

UR AD ( RuOF 

.END  WALL  1 

>  = 

6 . 

03 

OR AD ( ROOF 

.END  WALL  2 

)  = 

6. 

03 

OR AD ( ROOF 

?  F  LOOR 

)  = 

676 . 

cr 

r_l 

OR AD <  2 

IDE 

WALL 

1 

.SIDE  WALL  2 

)  = 

00 

OR AD  ( 2 

IDE 

WALL 

1 

.END  WALL  1 

)  = 

cr  cr 

ORAD  ( 2 

IDE 

WALL 

1 

.END  WALL  2 

)  - 

55 

OR  AD  ( 2 

IDE 

WALL 

1 

? FLOOR 

= 

216. 

39 

ORAEi  (2 

IDE 

WALL 

.END  WALL  1 

)  = 

55 

OR  AD  ( 2 

IDE 

WALL 

.END  WALL  2 

)  = 

55 

ORAD  ( 2 

IDE 

WALL 

4J 

. FLOOR 

}  = 

216. 

39 

ORAD ( END 

WALL  1 

.END  WALL  2 

)  = 

00 

ORAD  <  END 

WALL  1 

.  FLOOF; 

)  = 

126. 

53 

ORAD ( END 

WALL  2 

. F LuuR 

)  — 

1 26 . 

tr  z* 

TOTAL  HEAT  LOSS  BY  CONVECTION  AND  RADIATION  IS:  4174. 

MEAN  THERMAL  INSULATION*  .013176 
+STOP+ 


( K  /  W ) 
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HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  21 


ADDED  INSULATION  VALUES  (K/W) 
R  < ROOF  )—  . 035290 

R (SIDE  WALL  1  )=  .069430 

R(SIDE  WALL  2  >=  .069430 


R ( END  WALL  I  )=  .111450 
R ( END  WALL  2  )  =  .111450 
R (FLOOR  )=  .000000 


CONDITION  NUMBER 


2.  449 


INTERNAL  WALL  SURFACE  TEMPERATURES  (C) 


F (ROOF 

)  = 

-12. 26 

T  (SIDE  WALL  1 

)  = 

-14.19 

T (SIDE  WALL  2 

)  =£ 

-14.19 

T (END  WALL  1 

)  = 

-13. 99 

T (END  WALL  2 

)  = 

- 1 3 . 99 

T ( FLOOR 

)  = 

-39.41 

HEAT  LOSS  THROUGH  THE  WALLS  (W) 
i  RuuF 

Q ( S IDE  WALL  1 
Q ( S IDE  WALL  2 
Q (END  WALL  1 
Q (END  WALL  2 
Q ( FLOOR 

RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  <W> 


UR  AD  (  RuOF 

.SIDE  WALL  1 

)  = 

16.64 

UR'AD  \  RuuF 

.SIDE  WALL  2 

)  = 

16.  64 

ORAD ( ROOF 

-END  WALL  1 

)  = 

8. 64 

OR AD (ROOF 

-END  WALL  2 

)  - 

8.  64 

ORAD ( ROOF 

, FLOOR 

)  = 

692. 45 

ORAD ( S I DE 

WALL 

1 

. SIDE  WALL  2 

)  = 

.  00 

ORAD  ( 9.  [  DE 

WALL 

i 

-  END  WALL  1 

)  = 

-.54 

ORAD ( S I DE 

WALL 

1 

-  END  WALL  2 

>  = 

- .  54 

ORAD 1  SIDE 

WALL 

1 

•  FLOOR 

>  = 

217. 76 

□PAD  '’SIDE 

WALL 

.END  WALL  1 

>  = 

-.  54 

ORAD ■ SIDE 

WALL 

-  END  WALL  2 

>  - 

-.54 

OP AD ' S I DE 

WALL 

-  FLOOR 

)  = 

217. 76 

'.'RAD  (END1 

WALL  1 

-  END  WALL  2 

.  = 

.  00 

ORAD 'END 

WALL  1 

-  FLOOR 

)  = 

127.02 

ORAD  END 

WALL  2 

-  FLOOR 

i  = 

127. 02 

TOTAL  HEAT 

LOSS 

BY 

C  ON', ‘EC  T I  ON 

AND 

RADIATION  I 

)=  780.10 
)=  368.98 
)  =  368.98 
)=  231.10 
) =  23 1.10 
>=  2234.80 


MEAN  THERMAL  (NSULATION= 
*STOP+ 


49 


:>6  (W) 
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HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  EC- 


ADDED  INSULATION  VALUES  (K/W) 

R (ROOF  i=  .035290 

R ( S I DE  WALL  I  )=  .069430 

R ( S IDE  WALL  2  >=  .069430 

F:  (END  WALL  1  )=  .111450 

R  <  END  WALL  2  )=  .111450 

R (FLOOR  )=  .000000 

CONDITION  NUMBER  =  3.491 

INTERNAL  WALL  SURFACE  TEMPERATURES  <C> 


T ( ROOF 

)  = 

-IS. 16 

T (SIDE  WALL  1 

)  ~ 

-14. 92 

T (SIDE  WALL  2 

)  = 

-14. 92 

T(END  WALL  1 

)  = 

-14. 72 

T (END  WALL  2 

)  = 

-14.72 

T  FLOOR 

)  = 

-39 . 56 

HEAT  LOSS  THROUGH 

THE  WALL 

0 ( ROOF 

)  — 

614.22 

a (SIDE  WALL  1 

)  = 

358 . 59 

0 (SIDE  WALL  2 

)  = 

358. 59 

0 (END  WALL  1 

)  — 

224 . 60 

Q  (END  WALL  2 

)  = 

224 . 60 

0 ( FuOOR 

)  = 

1659. 74 

RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  (W) 


GRAD (ROOF 

, SIDE  WALL  1 

)  = 

-28. 01 

OR  A  Ci  ( RuOF 

..  SIDE  WALL  2 

»  = 

-28. 01 

OR AD (ROOF 

.END  WALL  1 

)  = 

-17. 18 

OR AD ( ROOF 

, END  WALL  2 

)  = 

-17. 18 

OR AD ( RUOF 

..  FLOOR 

)  = 

545.89 

OR AD ( S IDE  WALL  1 

..SIDE  WALL  2 

)  = 

.  00 

GRAD (SIDE  WALL  1 

..END  WALL  1 

)  = 

cr 

.  --0 

ORAD 'SIDE  WALL  1 

. END  WALL  2 

>  = 

53 

GRAD (SIDE  WALL  1 

, FLOOR 

)  = 

212. 30 

ORAD (SIDE  WALL  2 

?  END  WALL  1 

)  = 

- .  53 

ORAD (SIDE  WALL  2 

?  END  WALL  2 

)  = 

—  5  3 

ORAD (SIDE  WALL  2 

..  FLOOR 

)  = 

212. 80 

ORAD (END  WALL  1 

..END  WALL  2 

>  = 

.  00 

ORAD (END  WALL  1 

, FLOOR 

)  - 

124.12 

ORAD ( END  WALL  2 

. FLOOR 

\  = 

124.  12 

TOTAL  HEAT  LOSS  BY  CONVECTION  AND  RADIATION  IS:  344C 

MEAN  THERMAL  INSULATION=  . 015937 
♦STOP* 


( k  /'  w . 


HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC 
TEST  NUMBER :  1 1 


ADDED  INSULATION  VALUES  (K/W) 


R  ( R  u  U  F 

)  = 

.  000000 

R (SIDE  WALL  1 

)  = 

.  OOUOOO 

R < S IDE  WALL  2 

)  = 

. 000000 

R ( END  WALL  1 

)  = 

•  000000 

R ( END  WALL  2 

)  = 

. 000000 

R (FLOOR 

)  = 

•  000000 

CONDITION  NUMBER  =  1.801 

INTERNAL  WALL  SURFACE  TEMPERATURES  (C) 


T  (ROOF  :• 

-S3 . 90 

T (SIDE  WALL  1  ) 

= 

—  *  v .  *1  y 

T (SIDE  WALL  2  ) 

- 

— y  V  •  L 

T (END  WALL  1  ) 

= 

-39. 07 

T (END  WALL  2  ) 

=r 

-39. 07 

T (FLOOR  ) 

rr 

—39. 76 

HEAT  LOSS  THROUGH 

THE  WALL 

Q ( ROUF 

)  ~ 

4148. 06 

Q (SIDE  WALL  1 

)  -- 

1385. 06 

□(SIDE  WALL  2 

)  = 

1385.06 

□(END  WALL  1 

)  = 

850. 51 

□(END  WALL  2 

)  = 

850. 51 

Q ( FLOOR 

)  * 

895 . 53 

RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  (W) 


□RAD (ROOF 

S I  DE  WALL  l 

)  = 

3 . 30 

□RAD ( RUUF 

, SIDE  WALL  2 

)  = 

3 .  30 

□RAD (ROOF 

..END  WALL  1 

)  = 

86 

□RAD ( ROOF 

>  END  WALL  2 

)  = 

.  86 

□RAD ( ROOF 

..  FLOOR 

)  = 

21.99 

□RAD ( S I DE 

WALL 

1 

..SIDE  WALL  2 

)  * 

.  00 

□RAD ( S I DE 

WALL 

1 

, END  WALL  1 

)  = 

-.57 

□RAD (SIDE 

WALL 

1 

,END  WALL  2 

)  = 

- .  57 

□RAD ( S I DE 

WALL 

1 

, FLOOR 

)  = 

4.15 

□RAD < S I DE 

WALL 

, END  WALL  1 

)  = 

-.57 

□RAD  ( S I DE 

WALL 

zz 

?  END  WALL  2 

>  = 

- .  57 

□RAD ( S I DE 

WALL 

jZ 

, FLOOR 

)  = 

4.  15 

□RAD (END 

WALL  1 

..END  WALL  2 

)  = 

.  00 

□RAD (END 

WALL  1 

, FLOOR 

)  = 

3.  45 

□RAD ( END 

WALL  2 

.  FLOOR 

)  = 

3.45 

fOTAL  HEAT 

LOSS 

BY 

CONVECTION 

AND 

RADIATION 

MEAN  THERMAL  INSULATION*  .005781  (K/W) 

+-ST  0P+: 


HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  10 


ADDED  INSULATI 

JN 

VALUES  (K/W 

R  ( R  i„f  0  P 

)  = 

.  000000 

R (S IDE  WALL  1 

)  = 

. 000000 

R (SIDE  WALL  2 

)  = 

. 000000 

R ( END  WALL  1 

)  = 

. 000000 

R (END  WALL  2 

)  =r 

. 000000 

R (  FLUOR 

)  = 

. 000000 

CONDITION  NUMBER  =  2.145 

INTERNAL  WALL  SURFACE  TEMPERATURES  (C) 


T ( ROOF 

>  = 

-39. 53 

T  (SIDE  WALL  1 

)  = 

-39.44 

T (SIDE  WALL  2 

)  = 

-39. 44 

T ( END  WALL  1 

)  = 

-39 . 28 

T (END  WALL  2 

)  - 

-39. 28 

T ( FLOOR 

)  = 

-39 . 88 

HEAT  LOSS  THROUGH 

THE  WALL 

Q  ( ROOF 

)  = 

1754. 75 

Q ( S IDE  WALL  1 

)  = 

1 085 . 02 

Q(SIDE  WALL  2 

)  = 

1085. 02 

Q ( END  WALL  1 

)  = 

664 . 88 

Q  (END  WALL  2 

)  = 

664 . 88 

Q (FLOOR 

)  — 

464.87 

RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  (W) 


OR AD ( ROOF 

.-SIDE  WALL  1 

)  = 

-.  84 

OR AD (ROOF 

?  SIDE  WALL  2 

)  = 

-.84 

QRAD ( ROOF 

..END  WALL  1 

)  = 

- 1 . 30 

OR AD ( ROOF 

?  ENEi  WALL  2 

)  = 

- 1 . 30 

QRAD ( ROOF 

..  FLOOR 

)  = 

8.  72 

QRAD  ( 

SIDE 

WALL 

* 

X 

?  SIDE  WALL  2 

)  = 

.  00 

QRAD  ( 

SIDE 

WALL 

1 

..END  WALL  1 

>  = 

-.  44 

QRAD  ( 

S I DE 

WALL 

1 

?  END  WALL  2 

)  = 

-.44 

QRAD  ( 

3 1  DE 

WALL 

1 

..  FLOOR 

>  = 

3 . 79 

QRAD  ( 

SIDE 

WALL 

^2 

7  END  WALL  1 

)  = 

44 

QRAD  ( 

SIDE 

WALL 

..END  WALL  2 

)  = 

-.  44 

QRAD  ( 

5 1  DE 

WALL 

^2* 

..  FLOOR 

)  = 

3.  79 

QRAD (END  WALL  1 

..END  WALL  2 

)  = 

.  00 

QRAD (END  WALL  1 

. FLUOR 

)  = 

3.01 

QRAD (END  WALL  2 

, FLOOR 

)  = 

3 . 0 1 

TOTAL 

HEAT 

LOSS 

BY 

CONVECTION 

AND 

RADIATION  IS 

MEAN  THERMAL  INSULATION* 
♦STOP* 


009616  (K/W) 


45 


APPENDIX  D 


DATA  FROM  THEORETICAL  ANALYSIS 


44 


4  76  C 

4  77  CALL  SOLVE(NDIM,N,LU,NPIV,R,Z) 

4  78  C 

4  79  C  THE  IMPROVED  APPROXIMATION  IS  COMPUTED 

4  80  C 

4  81  DO  30  1=1, N 

4  82  X(I)=X(I)+Z(I) 

4  83-30  CONTINUE 

4  84  C 

4  85  RETURN 

4  86  END 


5  1C 
5  2  C 

5  3  FUNCTION  DTPROD(N,X,Y) 

5  4  C 

5  5  C 

5  6  C  THE  PURPOSE  OF  THIS  SUBPROGRAM  IS  TO  CALCULATE  THE  DOT 

5  7  C  PRODUCT  OF  TWO  VECTORS  OF  LENGTH  N.  THE  DOT  PRODUCT  IS 

5  8  C  ACCUMULATED  IN  DOUBLE  PRECISION  TO  ENSURE  MAXIMUM  ACCURACY. 

5  9  C 

5  10  C  CALLING  SEQUENCE:  SCALAR  CONSTANT  =  DTPROD(N,X,Y) 

5  11  C 

5  12  C  PARAMETERS: 

5  13  C 

5  14  C  N  -  AN  INTEGER  THAT  INDICATES  THE  LENGTH  OF  THE 

5  15  C  TWO  VECTORS 

5  16  C 

5  17  C  X,Y  -  REAL  VECTORS  OF  LENGTH  N 

5  18  C 

5  19  C 

5  20  C 

5  21  DOUBLE  PRECISION  XX,YY,SUM 

5  22  DIMENSION  X(50),Y(50) 

5  23  SUM  =  O.DO 

5  24  DO  10  1=1, N 

5  25  XX  =  X( I) 

5  26  YY  =  Y(  I) 

5  27  SUM  =  SUM  +  XX*YY 

5  28  10  CONTINUE 

5  29  DTPROD  =  SUM 

5  30  RETURN 

5  31  END 


IN  THE  CALLING  PROGRAM 

AN  INTEGER  CONSTANT  INDICATING  THE  NUMBER  OF 
UNKNOWNS  IN  THE  SYSTEM 

A  REAL  2-D  ARRAY  OF  SIZE  NDIM*N  HOLDING  THE 
ORIGINAL  MATRIX.  THIS  ARRAY  IS  NOT  ALTERED  BY 
IMPRV. 

A  REAL  2-D  ARRAY  OF  SIZE  NDIM*N  HOLDING  THE 
LU  DECOMPOSITON  OF  THE  ORIGINAL  MATRIX.  THIS 
ARRAY  IS  NOT  ALTERED  BY  IMPRV 

A  INTEGER  VECTOR  OF  SIZE  N  HOLDING  THE  PIVOT 
INFORMATION  FROM  THE  ELIMINATION  STEP. 

A  REAL  VECTOR  OF  SIZE  N  HOLDING  RIGHT  HAND 
SIDE  OF  THE  ORIGINAL  SYSTEM  TO  BE  SOLVED.  THIS 
ARRAY  IS  NOT  ALTERED  BY  IMPRV. 

A  REAL  VECTOR  OF  SIZE  N  HOLDING  THE  INITAL 
APPROXIMATE  SOLUTION.  ON  RETURN  THIS  VECTOR 
WILL  CONTAIN  THE  IMPROVED  APPROXIMATION 


4  50 

4  51 


A  REAL  VECTOR  OF  SIZE  N  WHICH  IS  DECLARED  BUT 
NOT  INITIALIZED  BY  THE  CALLING  PROGRAM.  ON 
RETURN  Z  WILL  CONTAIN  THE  CORRECTIONS  TO  THE 
GIVEN  APPROXIMATE  SOLUTION 

A  REAL  VECTOR  OF  SIZE  N  WHICH  IS  DECLARED  BUT 
NOT  INITIALIZED  BY  THE  CALLING  PROGRAM.  ON 
RETURN  "R"  WILL  CONTAIN  THE  RESIDUAL  R=B-AX 
WHERE  X  IS  THE  INITIAL  APPROXIMATE  SOLUTION. 
DOUBLE  PRECISION  ARITHMETIC  IS  USED  IN 
COMPUTING  "R" 


4  59  REAL*8  AA,XX,SUM 

4  60  DIMENSION  A(50,50) ,B(50) ,X(50) ,Z(50) ,R(50) ,NPIV(50) 

4  61  REAL  LU(50,50) 

4  62  C 

4  63  C  CALCULATE  THE  RESIDUALS 

4  64  C 

4  65  DO  20  1=1, N 

4  66  SUM  =  B( I) 

4  67  DO  10  J-l.N 

4  68  AA=A( I , J) 

4  69  XX=X(J) 

4  70  SUM=SUM-AA*XX 

4  71  10  CONTINUE 

4  72  R( I)  =  SUM 

4  73  20  CONTINUE 

4  74  C 

4  75  C  THE  RESIDUAL  SYSTEM  IS  SOLVED 


•-  ■v*  C  .  i  ■  •  • ' .  •  v*.  * r  - 


55  X ( 1 ) =B (KP I VOT ) 

56  DO  30  K=2,N 

57  KPIVOT  =  NPIV(K) 

58  KM1  =  K-l 

59  SUM  =  B(KPIVOT) 

60  DO  20  J=1,KM1 

61  SUM  =  SUM  -  LU(KPIVOT, J)*X( J) 

62  -20  CONTINUE 

63  X(K)  =  SUM 

64  30  CONTINUE 

65  C 

66  C  BACK  SUBSTITUTION  BEGINS 

67  C 

68  X ( N )  =  X(N)/LU(KP I VOT , N ) 

69  K  =  N 

70  DO  50  1=2, N 

71  KP1  =  K 

72  K  =  K-l 

73  KPIVOT  =  NPIV(K) 

74  SUM  =  X(K) 

75  DO  40  J  =  KP1  ,N 

76  SUM  =  SUM  -  LU(KPIVOT, J)*X( J) 

77  40  CONTINUE 

78  X(K)  =  SUM/LU(KPIVOT,K) 

79  50  CONTINUE 

80  C 

81  RETURN 

82  END 


4  2  C 

4  3  SUBROUTINE  IMPRV(NDIM,N,A,LU,NPIV,B,X,Z,R) 

4  4  C 

4  5  C 

4  6  C  GIVEN  AN  APPROXIMATE  SOLUTION  FOR  A  LINEAR  SYSTEM  OF  EQUATIONS 

4  1  C  THIS  SUBROUTINE  CARRIES  OUT  ONE  ITERATION  OF  THE  ITERATIVE 

4  8  C  IMPROVEMENT  PROCESS  FOR  COMPUTING  A  BETTER  APPROXIMATE  SOLUTION 

4  9  C  IT  IS  ASSUMED  THAT  BOTH  THE  ORIGINAL  MATRIX  AND  ITS  LU 

4  10  C  DECOMPOSITION  ARE  AVAILABLE.  THIS  SUBROUTINE  CAN  BE  USED 

4  11  C  IN  CONJUNCTION  WITH  THE  SUBROUTINES  DCOMP,  DCOMS  AND  SOLVE 

4  12  C  TO  FIND  THE  SOLUTION  OF  A  LINEAR  SYSTEM  OF  EQUATIONS. 

4  13  C 

4  14  C  TAKEN  FROM:  STUNT  MANUAL,  DEPT  OF  COMPUTER  SCIENCE,  U  OF  T 

4  15  C 

4  16  C  CALLING  SEQUENCE:  CALL  IMPRV(NDIM,N,A,LU,NPIV,B,X,Z,R) 

4  17  C 

4  18  C  PARAMETERS: 

4  19  C 

4  20  C 


NDIM 


AN  INTEGER  CONSTANT  INDICATING  THE  NUMBER  OF 
ROWS  IN  THE  ARRAYS,  "A"  AND  "LU”  AS  DECLARED 


1 

2  SUBROUTINE  SOLVE(NDIM,N,LU,NPIV,B,X) 

3  C 

4  C 

5  C  THIS  SUBROUTINE  PERFORMS  THE  FORWARD  AND  BACKWARD  SUBSTITUTIONS 

6  C  STEPS  IN  THE  SOLUTION  OF  A  SYSTEM  OF  LINEAR  EQUATIONS  AX=B. 

7  C  IT  ASSUMES  THAT  THE  TRIANGULAR  (OR  LU)  FACTORIZATION  OF  A  HAS 

8  C  ALREADY  BEEN  COMPUTED  BY,  SAY,  THE  ROUTINE  DCOMP  OR  DCOMS.  IF 

9  C  EITHER  ROUTINE  INDICATES  THAT  "A"  IS  SINGULAR,  THEN  THE  USE  OF 

10  C  SOLVE  MAY  PRODUCE  AN  OVERFLOW  INDICATION. 

11  C 

12  C 

13  C  REPRODUCED  FROM:  STUNT  MANUAL,  DEPT  OF  COMPUTER  SCIENCE,  U  OF  T 

14  C 

15  C 

16  C  CALLING  SEQUENCE:  CALL  SOLVE(NDIM,N,LU,NPIV,B,X) 

17  C 

18  C  PARAMETERS: 

19  C 

20  C  NDIM  -  AN  INTEGER  INDICATING  THE  NUMBER  OF  ROWS  IN  THE 

21  C  ARRAY  "A"  AS  DECLARED  IN  THE  CALLING  PROGRAM 

22  C 

23  C  N  -  AN  INTEGER  CONSTANT  INDICATING  THE  SIZE  OF  THE 

24  C  SYSTEM  TO  BE  SOLVED. 

25  C 

26  C  LU  -  A  REAL  2-D  ARRAY  OF  SIZE  NDIM*N  CONTAINING  THE 

27  C  LU  DECOMPOSITON  OF  A.  THIS  ARRAY  IS  NOT  ALTERED 

28  C  BY  SOLVE. 

29  C 

30  C  NPIV  -  AN  INTEGER  VECTOR,  OF  DIMENSION  N,  HOLDING  THE 

31  C  PIVOT  INFORMATION  FOR  THE  ELIMINATION  STEP. 

32  C 

33  C  X  -  A  REAL  VECTOR,  OF  SIZE  N,  THAT  IS  DECLARED  BUT 

34  C  NOT  INITIALIZED  BY  THE  CALLING  PROGRAM.  ON 

35  C  RETURN,  THIS  ARRAY  CONTAINS  THE  COMPUTED 

36  C  SOLUTION  OF  THE  SYSTEM. 

37  C 

38  C  B  -  A  REAL  VECTOR,  OF  SIZE  N,  HOLDING  THE  RIGHT  HAND 

39  C  SIDE  OF  THE  ORIGINAL  SYSTEM  TO  BE  SOLVED.  THE 

40  C  CONTENTS  OF  THIS  VECTOR  ARE  UNALTERED  BY  SOLVE. 

41  C 

42  C 

43  DIMENSION  B ( 50 ) ,X ( 50 ) , NP I V ( 50 ) 

44  REAL  LU(50,50) 

45  C 

46  C  CHECK  FOR  SYSTEM  OF  ONLY  ONE  UNKNOWN 

47  C 

48  IF  (N.GT.l)  GOTO  10 

49  X(l)  *  B(1)/LU{1 ,1) 

50  RETURN 

51  C 

52  C  FORWARD  ELIMINATION  ON  "B".  THE  RESULT  IS  PLACED  IN  X 

53  C 

54  10  KPIVOT  =  NPIV(l) 


40 


2 

91 

C 

2 

92 

DO  30  J=I ,N 

2 

93 

IP  *  NPIV(J) 

2 

94 

HOLD  =  ABS(A(IP,I))/D(IP) 

2 

95 

IF  ( HOLD. LE. COLMAX)  GOTO  30 

2 

96 

COLMAX  =  HOLD 

2 

97 

NROW  =  J 

2 

98 

■30 

CONTINUE 

2 

99 

C 

2 

100 

C 

TEST  FOR  SINGULARITY.  THE  MATRIX  IS  ASSUMED  TO  BE  SINGULAR 

2 

101 

C 

IF 

COLMAX  (THE  ABS.  VALUE  OF  THE  SCALED  PIVOT)  IS  EQUIVALENT 

2 

102 

C 

TO  ZERO,  IE,  1.0  +  COLMAX  =  1.0 

2 

103 

C 

2 

104 

c 

IF 

THIS  IS  TRUE  THEN  THE  ROUTINE  PROCEEDS  ON  TO  THE  (I+1)TH 

2 

105 

c 

STAGE  OF  THE  ELIMINATION 

2 

106 

c 

2 

107 

IF  ((l.O+COLMAX)  .NE.  1.0)  GOTO  40 

2 

108 

IND  =  -1 

2 

109 

GOTO  80 

2 

110 

c 

2 

111 

c 

IF 

AN  INTERCHANGE  IS  NECESSARY,  ALTER  THE  PIVOT  VECTOR  "NPIV" 

2 

112 

c 

2 

113 

40 

IPIVOT  =  NPIV(NROW) 

2 

114 

IF  (NROW  .EQ.  I)  GOTO  50 

2 

115 

NPIV (NROW)  »  NPIV(I) 

2 

116 

NPIV(I)  =  IPIVOT 

2 

117 

C 

2 

118 

C 

THE  MULTIPLIERS  FOR  THE  COMPUTATION  OF  THE  REMAINING  ROWS  ARE 

2 

119 

C 

DETERMINED  AND  ELINIMATION  IS  PERFORMED.  THE  VALUE  OF  EACH 

2 

120 

c 

MULTIPLIER  IS  STORED  IN  THE  POSITION  OF  THE  ELIMINATED 

2 

121 

c 

ELEMENT. 

2 

122 

c 

2 

123 

50 

IP1  =  1+1 

2 

124 

DO  70  J-IPl.N 

2 

125 

JPIVOT  =  NPIV(J) 

2 

126 

AMULT  =  A( JPIVOT, I) /A( IPIVOT, I) 

2 

127 

A( JPIVOT, I )  =  AMULT 

2 

128 

C 

2 

129 

DO  60  K=IP1,N 

2 

130 

A( JPIVOT, K)=A( JPIVOT, K)-AMULT*A( IPIVOT, K) 

2 

131 

60 

CONTINUE 

2 

132 

70 

CONTINUE 

2  133  80  CONTINUE 

2  134  C 

2  135  TEST  -  1.  +  ABS(A(NPIV{N),N)) 

2  136  IF  (TEST  .EQ.  1.)  THEN 

2  137  IND  =  -1 

2  138  END  IF 

2  139  C 

2  140  RETURN 

2  141  END 


A 


A  REAL  2-DIMENSIONAL  ARRAY,  OF  SIZE  NDIM*N 
HOLOING  THE  MATRIX  TO  BE  DECOMPOSED.  ON 
RETURN,  THE  CONTENTS  OF  "A"  ARE  REPLACED 
BY  THE  LU  FACTORIZATION. 

AN  INTEGER  VECTOR,  OF  SIZE  N,  WHICH  IS 
UNINITIALIZED  AT  THE  TIME  OF  CALLING.  THIS 
ARRAY  WILL  RECORD  THE  REARRANGING  OF  THE 

A  REAL  VECTOR,  OF  SIZE  N,  THAT  IS  USED  AS 
A  WORKSPACE  FOR  THE  SCALING  OPERATION.  THIS 
ARRAY  IS  DECLARED  BUT  NOT  INITIALIZED  IN 
THE  CALLING  PROGRAM. 

AN  INTEGER  INDICATING  IF  "A"  IS  SINGULAR  OR  NOT 
IND  =  -1  ;  "SCALED"  PIVOT  >  UNIT  ROUND-OFF 
=  0  ;  "A"  IS  NONSINGULAR 

THIS  PARAMETER  IS  NOT  INITIALIZED  AT  THE  TIME 
OF  CALLING. 


DIMENSION  A(50,50) ,0(50) ,NPIV(50) 

IND  =  0 

CHECH  FOR  A  SYSTEM  OF  ONLY  ONE  UNKNOWN 

IF  (N.EQ.l)  RETURN 

INITIALIZE  PIVOT  AND  0  VECTORS 

DO  20  1=1, N 
NPIV(I)  =  I 
D( I ) =0 . EO 

THE  LARGEST  ABSOLUTE  VALUES  IN  EACH  ROW  ARE  RECORDED  IN  "D" 
DO  10  J=1,N 

D( I)  =  AMAX1(D( I) ,ABS(A(  I , J) ) ) 

CONTINUE 

IF  (D(I).EQ.O.EO)  D( I)  =  l.EO 
CONTINUE 

MAIN  LOOP  FOR  GAUSS  ELIMINATION 
NM1  =  N-l 
DO  80  I=1,NM1 

DETERMINE  THE  LARGEST  "SCALED  PlV^r  ie, 

MAX  A( J, I)/D( J)  ,  I»J>-N 


S 


38 


1  217  270  CONTINUE 

1  218  C 

1  219  280  FORMAT ( '  QRAD( ' ,A12, ' , 1 ,A12, ' )=  ',F8.2) 

1  220  C 

1  221  WRITE (*,290)  QTOT 

1  222  290  FORMAT(/, '  TOTAL  HEAT  LOSS  BY  CONVECTION  AND  RADIATION  ' 

1  223  +' IS:  ' ,F10.2, '  (W) ‘ ) 

1  224  C 

1  225  WRITE (*,300)  (TIN  -  T0UT)/QT0T 

1  226  300  F0RMAT(/, '  MEAN  THERMAL  INSULATION=  ',F10.6,'  (K/W) ' ) 

1  227  C 

1  228  RETURN 

1  229  END 


2  1  C 

2  2  C 

2  3  SUBROUTINE  DCOMS(NDIM,N,A,NPIV,D, IND) 

2  4  C 

2  5  C 

2  6  C  THIS  SUBROUTINE,  DECOMPOSITION  WITH  SCALED  PARTIAL  PIVOTING, 

2  7  C  DOES  GAUSSIAN  ELIMINATION  OR,  EQUIVALENTLY,  A  TRIANGULAR  (LU) 

2  8  C  FACTORIZATION  OF  THE  N*N  MATRIX  STORED  IN  THE  ARRAY  "A".  AT 

2  9  C  COMPLETETION,  THE  "A"  WILL  CONTAIN  THE  LOWER  TRIANGULAR  MATRIX 

2  10  C  OF  MULTIPLIERS  USED  IN  THE  ELIMINATION  AS  WELL  AS  THE  UPPER 

2  11  C  TRIANGULAR  MATRIX  "U",  THE  RESULT  OF  THE  ELIMINATION. 

2  12  C 

2  13  C  THE  MATRIX  IS  ASSUMED  TO  BE  SINGULAR  IF  EITHER  SOME  ROW  IS 

2  14  C  ZERO  INITIALLY,  OR,  SOME  "SCALED"  PIVOT  DURING  THE  ELIMINATION 

2  15  C  IS  SMALLER  THAN  UNIT  ROUND-OFF.  IF  THE  FORMER  HOLDS,  THE 

2  16  C  DECOMPOSITION  DOES  NOT  COMMENCE.  IN  THE  LATTER  CASE,  DCOMS 

2  17  C  WILL  COMPLETE  THE  DECOMPOSITION  BUT  THE  RESULTING  UPPER 

2  18  C  TRIANGULAR  MATRIX  WILL  BE  SINGULAR. 

2  19  C 

2  20  C  THIS  ROUTINE  CAN  BE  USED  IN  CONJUNCTION  WITH  THE  ROUTINE 
2  21  C  SOLVE  TO  FIND  THE  SOLUTION  TO  A  SYSTEM  OF  LINEAR  EQUATIONS. 

2  22  C  THE  CURRENT  MAXIMUM  SIZE  OF  THE  SYSTEM  IS  50*50 

2  23  C 

2  24  C  TAKEN  FROM:  STUNT  MANUAL;  DEPT  OF  COMPUTER  SCIENCE;U  OF  T 

2  25  C 

2  26  C 

2  27  C  CALLING  SEQUENCE:  CALL  DCOMS(NDIM,N,A,NPIV,D, IND) 

2  28  C 

2  29  C  PARAMETERS: 

2  30  C 

2  31  C  NDIM  -  AN  INTEGER  INDICATING  THE  NUMBER  OF  ROWS 

2  32  C  IN  THE  ARRAY  "A"  AS  DECLARED  IN  THE 

2  33  C  CALLING  PROGRAM. 

2  34  C 

2  35  C  N  -  AN  INTEGER  CONSTANT  INDICATING  THE  SIZE  OF 

2  36  C  THE  SYSTEM  TO  BE  SOLVED 


1  163  IF  (IND)  120,140,140 

1  164  120  WRITE(*,130) 

1  165  130  FORMAT ( '  ' , 'MATRIX  IS  SINGULAR,  TOO  BAD!  PROCESS  ABORTED') 

1  166  STOP 

1  167  C 

1  168  140  CALL  SOLVE(NDIM,N,LU,NPIV,B,T) 

1  169  CALL  IMPRV(NDIM,N,A,LU,NPIV,B,T,Z,R) 

1  170  C 

1  171  C  CONDA  GIVES  A  LOWER  BOUND  FOR  THE  CONDITION  NUMBER 

1  172  C  OF  THE  MATRIX 

1  173  C 

1  174  CONDA=  SQRT ( DTP  ROD ( N , Z , Z ) /DTP  ROD ( N , T , T ) ) *16 . **6 

1  175  WRITE(*,150)  CONDA 

1  176  150  FORMAT ( 1  * , • CONDITION  NUMBER  =  * ,F15.3) 

1  177  C 

1  178  WRITE(*,160) 

1  179  160  FORMAT(/, 1  INTERNAL  WALL  SURFACE  TEMPERATURES  (C) ' ,/) 

1  180  C 

1  181  DO  170  1=1, N 

1  182  WRITE(*,180)  AREA(I),T(I) 

1  183  170  CONTINUE 

1  184  180  FORMAT ( '  * ,'T( ' ,A12, ')*  ',F8.2) 

1  185  C 

1  186  C  NOW  CALCULATE  THE  HEAT  TRANSFER  RATES  THROUGH  THE  WALLS 

1  187  C  TO  THE  OUTSIDE  AIR,QCON(I),  AND  THE  TOTAL  HEAT  LOSS  RATE,  QTOT 

1  188  C 

1  189  QTOT  =0.0 

1  190  DO  200  1=1,6 

1  191  QCON(I)  =  (T( I)  -  TOUT)/RESIST( I) 

1  192  QTOT  =  QTOT  +  QCON(I) 

1  193  200  CONTINUE 

1  194  C 

1  195  C  NOW  CALCULATE  THE  RADIANT  HEAT  TRANSFER  BETWEEN  THE  WALLS 
1  196  C 

1  197  DO  220  1=1,6 

1  198  DO  210  J=I+1,6 

1  199  QRAD( I , J)  =  (T(I)-T(J))/RRAD(I,J) 

1  200  QRAD( J, I)  =  -QRAD(I, J) 

1  201  210  CONTINUE 

1  202  220  CONTINUE 
1  203  C 

1  204  WRITE ( *, 230 ) 

1  205  230  FORMAT(/, '  HEAT  LOSS  THROUGH  THE  WALLS  (W) ’ ,/) 

1  206  C 

1  207  WRITE (*,240)  ( ( AREA( I ) ,QC0N( I) ) , 1=1,6) 

1  208  240  FORMAT ( '  V  Q(',A12,')=  ’,F8.2) 

1  209  C 

1  210  WRITE(*,250) 

1  211  250  FORMAT{/, '  RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  (W) ’ ,/) 

1  212  C 

1  213  DO  270  1=1,6 

1  214  DO  260  J= 1+1,6 

1  215  WRITE(*,280)  AREA( I) ,AREA( J) ,QRAD( I, J) 

1  216  260  CONTINUE 


H  E  A  T  TRANSft  E  ft  C  A  L  U  U  L  A  f  1 U  N S  F  U  ft;  A  N  I  DEAL  i  ZEE1  A  ft'  C  MODEL 
TEST  NUMBER:  30 


ADDED  INSULATION  VALUES  <K/W> 
R ( ROOF  )=  .033300 

ft' (SIDE  WALL  1  )=  .069400 

ft;  ( SIDE  WALL  V  ) —  .069400 

R (END  WALL  1  )=  .111500 

ft:  (END  WALL  2  )=  .111500 

ft;  (FLOOR  ) =  . 052000 

CONDITION  NUMBER  = 


3.  1 58 


INTERNAL  WALL  SURFACE  TEMPERATURES  (C) 


T  ( R'UUft  —  —  1 1 . 38 
T (SIDE  WALL  1  )=  -9.69 
T (SI DE  WALL  2  >=  -9.69 
T  (END  WALL  1  >=  -9.  £.7 
T (END  WALL  2  )=  -9.67 
T (FLOOR  )-  -14.99 

HEAT  LOSS  THROUGH  THE  WALLS  (W) 

Q ( ROOF  )-  804.86 
Q ( S IDE  WALL  1  )=  433.54 
Q  ( S'  IDE  WALL  2  )=  433.54 
0(END  WALL  1  )=  269.38 
Q (END  WALL  2  )=  269.38 
GKFLOOF;  )=  478.59 


QRAD ( ROOF 
QRAD ( ROOF 
QRAD (ROOF 
QRAD (ROOF 
QRAD ( ROOF 
QRAD (SIDE  WALL  : 
QRAD (SIDE  WALL  : 
QRAD (SIDE  WALL  ; 
QRAD (SIDE  WALL  : 
QRAD (SIDE  WALL  ( 
QRAD (SIDE  WALL  ( 
QRAD (SIDE  WALL  ( 
QRAD (END  WALL  1 
QRAD (END  WALL  1 
QRAD ( END  WALL  2 


FER  BETWEEN 

WALLS 

(W) 

» SIDE  WALL  1 

)  = 

-14. 57 

.-SIDE  WALL  2 

)  = 

-14. 57 

E!  WALL  1 

>  = 

-8 . 52 

?  END  WALL  2 

)  = 

-8.  52 

FLOOR 

)  = 

92.  12 

-  SIDE  WALL  2 

)  = 

.  00 

?  END  WALL  1 

)  = 

-.  05 

.-END  WALL  2 

)  = 

-.05 

FLUOR 

)  = 

45 . 75 

.-END  WALL  1 

>  = 

- .  05 

.-END  WALL  2 

)  = 

-.  05 

- FLUUR 

)  = 

45.  75 

.-END  WALL  2 

)  = 

.  00 

FLOOR 

)  = 

26.  57 

-  FLUOR 

)  = 

26.  57 

TOTAL  HEAT  LOSS  BY  CONVECTION  AND  RADIATION  I: 

MEAN  THERMAL  INSULATION*  .020451  (K/Wi 
+STUP* 


Y*  rf*  I 


52 


a 


HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  31 


iDDED  INSULATION  VALUES  (K/W) 


E(ROOF  )~ 

R(SlDE  WALL  1  )— 

R (side  wall  l  )  — 

R (END  WALL  1  )= 

F'  (END  WALL  2  )  = 

R( FLOOR  ')- 

CONDITION  NUMBER  = 


. 035300 
. 069400 
. 069400 
.111 500 
. 1 1 1500 
. 052000 


4.050 


INTERNAL  WALL  SURFACE  TEMPERATURES  <C> 


r  (  r  u  u  i~ 

T( SIDE  WALL  1 
TCSIDE  WALL  2 
T (END  WALL  1 
T  (END  WALl.  2 
T \ FLOOR 


LOSS 


)  =  -11.97 

THROUGH  THE  WALLS  (W) 


Q  ( ROOF  )=  931.87 
Q  (SIDE  WALL  1  )=  440.23 
Q  < S I DE  WALL  2  )=  440.23 
Q  (END  WALL  1  )=  273.36 
Q(END  WALL  2  )=  273.36 


Q ( FLOOR 

) 

536. 30 

(ADI ANT  HEAT  TRANS 

•FER  BETWEEN 

WALLS 

(W) 

QRAD ( ROOF 

.-SIDE  WALL  1 

)  = 

20.  40 

QRAD ( ROOF 

, SIDE  WALL  2 

)  = 

20 .  40 

QRAD  ( RiJUF 

-  END  WALL  1 

)  = 

1 1 . 32 

QRAD ( ROOF 

-  END  WALL  2 

)  = 

1 1 . 32 

QRAD ( RUUF 

FLOOR 

)  = 

130. 41 

QRAD (SIDE 

WALL 

1 

, SIDE  WALL  2 

)  = 

.  00 

QRAD ( S I DE 

WALL 

1 

.-END  WALL  1 

)  = 

.01 

QRAD ( S I DE 

WALL 

1 

?  END  WALL  2 

)  = 

.  01 

QRAD (SIDE 

WALL 

1 

-  FLOOR 

)  = 

23.  75 

QRAD (SIDE 

WALL 

2 

?  END  WALL  1 

>  = 

.  01 

QRAD (SIDE 

WALL 

?  END  WALL  2 

)  = 

.  01 

QRAD (SIDE 

WALL 

, FLOOR 

)  = 

23.  75 

QRAD (END  WALL  1 

?  END  WALL  2 

)  = 

.  00 

QRAD (END  i 

WALL  1 

.-  FLOOR 

>  = 

13.73 

QRAD (END  WALL  2 

-  FLOOR 

)  = 

13.  73 

TOTAL  HEAT  LOSS  BY  CONVECTION  AND  RADIATION  IS: 

MEAN  THERMAL  INSULATION®  .018996  (K/W) 

♦STOP* 


■895 . 34  (W 


»  v  -•*  \ 


HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  .32 


ADDED  INSULATION  VALUES  (K/W) 


R  RijijR 

)  = 

. 033300 

R  (SIDE  WALL  1 

)  rr 

. 069400 

RESIDE  WALL  2 

)  = 

. 069400 

R  (END  WALL  1 

)  = 

.111 500 

R (END  WALL  2 

)  = 

.111 500 

R (FLOOR 

)  = 

. 052000 

CONDITION  NUMBER  = 

3.  S31 

INTERNAL  WALL 

SURFACE  TEMPERATURE: 

T (ROOF 

= 

—7.  39 

T (SIDE  WALL  1 

>  = 

-9 .  38 

T (SIDE  WALL  2 

)  = 

-9 . 38 

T  END  WALL  1 

)  = 

—9 . 38 

T (END  WALL  2 

>  = 

—9  _  33 

T  ( FLOOR! 

)  =r 

-12.36 

HEAT  LOSS  THROUGH 

THE  WALLS  (W) 

Q ( ROOF 

)  = 

9 1 6 . 36 

Q (SIDE  WALL  1 

)  = 

437. 96 

Q ( S IDE  WALL  2 

)  = 

437.96 

Q (END  WALL  1 

)  = 

272. 00 

Q  (END  WALL  2 

)  = 

272.00 

Q ( FLOOR 

)  = 

523. 92 

RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  (W) 


O.RAD  (ROOF 

, S I DE  WALL  1 

)  = 

17.16 

QRAD (ROOF 

.-SIDE  WALL  2 

)  = 

17.  16 

QRAD (ROOF 

-END  WALL  1 

)  = 

9.91 

QRAD ( ROOF 

.-END  WALL  2 

)  = 

9.  91 

QRAD ( ROOF 

.-  FLOOR 

)  = 

126.63 

QRAD (SIDE 

WALL  1 

.-  SIDE  WALL  2 

)  = 

.  00 

QRAD (SIDE 

WALL  1 

.-END  WALL  1 

)  = 

-.01 

QRAD (SIDE 

WALL  1 

.-END  WALL  2 

)  = 

-.01 

QRAD (SIDE 

WALL  1 

-  FLUOR.’ 

)  = 

25.  71 

QRAD (SIDE 

WALL  2 

-  END  WALL  1 

)  = 

- .  0 1 

QRAD ( S I DE 

WALL  2 

.-END  WALL  2 

)  = 

- .  0 1 

QRAD (SIDE 

WALL  2 

-  FLOOR 

)  = 

25.  71 

QRAD (END 

WALL  1 

.-END  WALL  2 

)  = 

.  00 

QRAD ( END 

WALL  1 

-  FLOOR 

)  = 

1 4 . 89 

QRAD ( END 

WALL  2 

.-  FLOOR; 

)  = 

1 4 . 89 

TOTAL  HEAT 

LOSS  BY 

CONVECTION 

AND 

RADIATION 

MEAN  THERMAL  INSULATION® 
+STGP* 


019193  (K/W) 


nt AT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  40 


ADDED  INSULATION  VALUES  (K/W) 

R (ROOF  > —  , 070600 

R (SIDE  WALL  1  )=  .138700 

R (SIDE  WALL  2  >=  .138700 

R ( END  WALL  1  )=  .222900 

R ( END  WALL  2  )=  .222900 

R (FLOOR  )=  .104000 

CONDITION  NUMBER  =  6.376 

INTERNAL  WALL  SURFACE  TEMPERATURES  (C) 

T ( ROOF  >=  -2.44 

T (SIDE  WALL  1  • =  -1.17 

T (SIDE  WALL  2  >  =  -1.17 

T'END  WALL  1  >  =  -1.15 

T  (END  WALL  2  .' =  -1.15 

T  (FLOOR  )=  -5.06 

HEAT  LOSS  THROUGH  THE  WALLS  (W) 

!.!  [  RlJUF  )  — 

ill  ( S  I DE  WALL  1  ;  = 

CM. SIDE  WALL  2  )  = 

Q (END  WALL  1 
0 (END  WALL  2  )= 

0 (FLOOR  }= 


RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  '(W) 


QRAD ( ROOF 

•SIDE  WALL  1 

)  = 

-10.93 

QRAD  ( 

ROOF 

•SIDE  WALL  2 

)  = 

- 10. 93 

QRAD ( ROOF 

•END  WALL  1 

)  = 

-6 . 40 

QRAD (ROOF 

•END  WALL  2 

)  = 

-6.40 

QRAD ( ROOF 

•  FLOOR 

)  = 

67 . 0 1 

QRAD  ( 

SIDE 

WALL 

1 

•SIDE  WALL  2 

)  = 

.  00 

QRAD  ( 

SIDE 

WALL 

1 

•END  WALL  1 

}  = 

-.04 

QRAD  ( 

SIDE 

WALL 

1 

•END  WALL  2 

)  = 

-.  04 

QRAD  0 

S I  DE 

WALL 

1 

•  FLOOR 

)  = 

33. 62 

QRAD  ( 

SIDE 

WALL 

£ 

•END  WALL  1 

)  = 

-.04 

QRAD  ( 

SIDE 

WALL 

•END  WALL  2 

>  = 

-.04 

QRAD  ( 

3 1  DE 

WALL 

•  FLOOR 

l  = 

QRAD (END  WALL  1 

•END  WALL  2 

>  = 

.  00 

QRAD (END  1 

WALL  1 

•  FLOOR 

>  = 

19.53 

QRAD (END  WALL  2 

•  FLOOR 

)  = 

19.53 

TOTAL 

HEAT 

LOSS 

BY 

CONVECTION 

AND 

RADIATION 

MEAN  THERMAL  INSULATION=  .031076  (K/W) 
♦STOP* 


530. 09 
278. 92 
278. 92 
1 73. 43 
173. 43 
335 .  09 


$£at  transfer  calculations  for  an  idealized  arc  model 

TEST  NUMBER:  41 


ADDED  INSULATION  VALUES  (K/W) 


R ( ROOF 

)  = 

. 070600 

R ( S IDE  WALL  1 

)  =r 

. 138700 

R ( SIDE  WALL  2 

)  = 

.  138700 

R  i.  END  WALL  1 

)  = 

. 222900 

R (END  WALL  2 

)  — 

. 222900 

R ( FLOOR 

)  — 

. 104000 

CONDITION  NUMBER  =  3.097 

INTERNAL  WALL  SURFACE  TEMPERATURES  (C) 


1  k  R  U  u  F 

)  = 

~7 

.  /  .*« 

T (SIDE  WALL  1 

)  = 

- 1 . 84 

T (SI DE  WALL  2 

)  = 

-1.84 

T (END  WALL  1 

)  = 

- 1 . 86 

T (END  WALL  2 

)  = 

- 1 . 86 

T (FLOOR 

)  = 

-3.41 

HEAT  LOSS  THROUGH  THE  WALLS  (W> 

Q  (ROOF  >  =  574. SI 

QCSIDE  WALL  1  )=  274.11 

Q (SIDE  WALL  2  )=  274.11 

Q (END  WALL  1  )=  170.26 

Q (END  WALL  2  )=  170.26 

Q (FLOOR  )=  350.95 

RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  (W) 


QRAD ( ROOF 

.SIDE  WALL  1 

)  = 

22.21 

QRAD (ROOF 

.  SIDE  WALL  2 

)  = 

22.  21 

QRAD (ROOF 

.END  WALL  1 

)  = 

12.  98 

QRAD (ROOF 

.END  WALL  2 

)  = 

1 2 .  98 

QRAD ( ROOF 

?  FLOOR 

)  = 

105. 66 

QRAD (SIDE 

WALL 

1 

.SIDE  WALL  2 

)  = 

.  00 

QRAD (SIDE 

WALL 

1 

.END  WALL  1 

>  = 

.  07 

QRAD (SIDE 

WALL 

1 

.END  WALL  2 

)  = 

.  07 

QRAD (SIDE 

WALL 

1 

.  FLOOR 

)  = 

1 3 . 56 

QRAD (SIDE 

WALL 

jL 

.END  WALL  1 

)  = 

.  07 

QRAD ( 3 I DE 

WALL 

.END  WALL  2 

i  = 

.  07 

QRAD (SIDE 

WALL 

. FLOOR 

)  = 

13.  56 

QRAD (END 

WALL  1 

.END  WALL  2 

)  = 

.  00 

QRAD ( END 

WALL  1 

?  FLOOR 

>  = 

7.  72 

QRAD (END 

WALi-  2 

. FLOOR 

)  = 

7.72 

TOTAL  HEAT 

LOSS 

BY 

CONVECTION 

AND 

RADIATION 

MEAN  THERMAL  INSULATION®  .030312  (K/W) 
♦STOP* 


56 


S' 


HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  42 


INSULATION  VALUES  (K/W) 


ADDED 
R ( ROOF 

R ( S I DE  WALL  1 
R ( S I DE  WALL  2 
R(END  WALL  1 
R (END  WALL  2 
R  <  FLOOR 


>  =  . 070600 

)  =  .138700 

>=  .138700 

>  =  . 222900 

>  --  .  222900 

)  =  .104000 


CONDITION  NUMBER 


3.  135 


INTERNAL  WALL  SURFACE  TEMPERATURES  <C> 
T ( ROOF 

T (SIDE  WALL  1 
T (SI DE  WALL  2 
T ( END  WALL  1 
T (END  WALL  2 
T (FLOOR 

HEAT  LOSS  THROUGH  THE  WALLS  (W) 


Q ( ROOF  ) = 

552 . 05 

Q ( S IDE  WALL  1  )= 

271. 40 

1 

Q ( S IDE  WALL  2  )= 

271.40 

« 

Q ( END  WALL  1  )= 

168.66 

Q  (END  WALL  2  )  = 

168.66 

Q (FLOOR  )= 

344 . 59 

i 

RADIANT  HEAT  TRANS 

FER  BETWEEN 

WALLS 

(W) 

Wi 

QRAD ( ROuF 

-SIDE  WALL  1 

)  = 

1 1 . 55 

QRAD ( ROOF 

-SIDE  WALL  2 

)  = 

1 1 . 55 

QRAD (ROOF 

-END  WALL  1 

)  = 

6.  71 

OR AD (ROOF 

-END  WALL  2 

)  = 

6.71 

QRAD (ROOF 

-  FLOOR 

)  = 

81.44 

ft 

QRAD (SIDE  WALL  1 

-SIDE  WALL  2 

)  = 

.  00 

QRAD (SIDE  WALL  1 

-END  WALL  1 

)  = 

.  01 

/■ 

QRAD (SIDE  WALL  1 

-END  WALL  2 

)  = 

.  01 

QRAD (SIDE  WALL  1 

-  FLOOR 

)  = 

16.  02 

QRAD (SIDE  WALL  2 

-END  WALL  1 

)  = 

.  01 

*■ 

QRAD (SIDE  WALL  2 

-END  WALL  2 

)  = 

.  0 1 

» 

QRAD (SIDE  WALL  2 

-  FLOOR 

)  = 

16.02 

QRAD (END  WALL  1 

-END  WALL  2 

)  = 

.  00 

QRAD (END  WALL  1 

?  FLOOR 

)  = 

9 .  25 

QRAD (END  WALL  2 

-  FLOOR 

}  = 

9 . 25 

TOTAL  HEAT  LOSS  BY  CONVECTION  AND  RADIATION  IS:  1776.77  (W) 

MEAN  THERMAL  INSULATION=  .030955  (K/W) 

•;  *sTop+ 


I 
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HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  APC  MODEL 
TEST  NUMBER:  50 


ADDED  INSULATION  VALUES  (K/W> 


R (ROOF 

)  = 

.141 200 

R (SIDE  WALL  1 

)  = 

. 277400 

R ( S I DE  WALL  2 

}  = 

. 277400 

R ( END  WALL  1 

)  = 

. 445800 

R ( END  WALL  2 

)  = 

.  445800 

R (FLOOR 

)  ~ 

. 207900 

CONDITION  NUMBER  =  1.220 


INTERNAL  WALL  SURFACE  TEMPERATURES  (C) 


T (ROOF  ) =  4.61 
T (SIDE  WALL  1  >=  5.42 
T (SIDE  WALL  2  )=  5.42 
T (END  WALL  1  )=  5.43 
T (END  WALL  2  )=  5.43 
T (FLOOR  >=  2.96 


HEAT  LOSS 

THROUGH 

THE  WALLS  <W) 

Q ( ROOF 

)  = 

315. 33 

Q ( S IDE  WALL  1  ) = 

163. 43 

Q ( S IDE  WALL  2  )= 

163. 43 

Q (END  WALL  1  )= 

101. 66 

O!  :  END  WALL  2  )  = 

101.66 

Q ( FLOOR 

)  = 

206. 37 

RADIANT  HEAT  TRAN 

5FER  BETWEEN,  WALLS 

(W) 

QRAD (ROOF 

-SIDE  WALL  1 

)  = 

-7.00 

QRAD ( ROOF 

.-SIDE  WALL  2 

)  = 

-7.  00 

OR AD ( ROOF 

.-END  WALL  1 

)  = 

-4.  10 

QRAD ( ROOF 

.-END  WALL  2 

)  = 

-4.  10 

QRAD ( ROOF 

.-  FLOOR 

)  = 

42 . 08 

QRAD (SIDE 

WALL  1 

.-SIDE  WALL  2 

)  = 

.  00 

QRAD (SIDE 

WALL  1 

-END  WALL  1 

)  = 

- .  03 

QRAD (SIDE 

WALL  1 

-END  WALL  2 

)  = 

- .  03 

QRAD ( S I DE 

WALL  1 

-  FLOOR 

)  = 

21 . 25 

QRAD (SIDE 

WALL  2 

-  END  WALL  1 

)  = 

- .  03 

QRAD (SIDE 

WALL  2 

-END  WALL  2 

)  = 

- .  03 

QRAD (SIDE 

WALL  2 

-  FLOOR 

)  = 

21.25 

QRAD ( END 

WALL  1 

-END  WALL  2 

)  = 

.  00 

QRAD (END 

WALL  1 

-  FLOOR 

)  = 

12.  35 

QRAD (END 

WALL  2 

-  FLOOR 

)  = 

12.  35 

TOTAL  HEAT  LOSS  BY  CONVECTION  AND  RADIATION 


052288  (K/W) 


1051 .  : 


i.-!. 


57  (W) 


MEAN  THERMAL  INSULATION3 
+STOP* 
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HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  51 


ADDED  INSULATION  VALUES  (K/W) 


R ( R u OF 

)  =5 

. 141200 

R  ( S IDE  WALL  1 

)  = 

.  277400 

R < S IDE  WALL  2 

)  = 

. 277400 

R  (END  WALL  1 

)  = 

. 445800 

R (END  WALL  2 

)  = 

.  445800 

R (FLOOR 

)  = 

. 207900 

CONDITION  NUMBER  =  2.449 

INTERNAL  WALL  SURFACE  TEMPERATURES  <C> 


T < ROOF 

)  = 

4 . 30 

T (SIDE  WALL  1 

)  = 

3.  19 

T(SIDE  WALL  2 

)  s 

3.  19 

T(END  WALL  1 

)  = 

3.  17 

T (END  WALL  2 

)  = 

3.  17 

T ( FLOOR 

)  = 

2.22 

HEAT  LOSS  THROUGH 

THE  WALL: 

Q (ROOF 

)  “ 

313. 15 

Q(SIDE  WALL  1 

)  = 

155. 40 

Q ( S IDE  WALL  2 

)  = 

155. 40 

Q (END  WALL  1 

)  ■= 

96.61 

Q(END  WALL  2 

)  = 

96.61 

Q ( FLOOR 

)  = 

202. 84 

RADIANT  HEAT  TRANSFER  BETWEEN  WALLS  (W> 


QRAD ( ROOF 

.SIDE  WALL  1 

>  = 

9.  60 

QRAD ( ROOF 

.SIDE  WALL  2 

)  = 

9.  60 

OR AD ( ROOF 

.END  WALL  1 

)  = 

5.62 

QRAD ( ROOF 

.END  WALL  2 

)  = 

5.  62 

QRAD (ROOF 

. FLOOR 

)  = 

52 . 98 

QRAD (SIDE 

WALL  1 

.SIDE  WALL  2 

)  = 

.  00 

QRAD (SIDE 

WALL  1 

.END  WALL  1 

)  = 

.04 

QRAD (SIDE 

WALL  1 

.END  WALL  2 

)  = 

.  04 

QRAD (SIDE 

WALL  1 

. FLOOR 

)  = 

8 .  33 

QRAD (SIDE 

WALL  2 

.END  WALL  1 

)  = 

.  04 

QRAD ( S I DE 

WALL  2 

.END  WALL  2 

)  = 

.  04 

QRAD (SIDE 

WALL  2 

?  FLUOR 

)  = 

8 . 33 

QRAD (END  1 

WALL  1 

.END  WALL  2 

)  = 

.  00 

QRAD (END  1 

WALL  1 

. FLOOR 

)  = 

4 . 76 

QRAD (END  1 

WALL  2 

. FLOOR 

)  = 

4 . 76 

TOTAL  HEAT 

LOSS  BY 

CONVECTION 

AND 

RADIATION 

MEAN  THERMAL  INSULATION^  .053921  (K/W) 
*STOP* 
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HEAT  TRANSFER  CALCULATIONS  FOR  AN  IDEALIZED  ARC  MODEL 
TEST  NUMBER:  52 


ADDED  INSULATION  VALUES  (K/W) 


R ( ROOF 

)  = 

.141 200 

R ( S I DE  WALL  1 

>  = 

. 277400 

R ( S I DE  WALL  2 

)  = 

.  277400 

R(END  WALL  1 

)  = 

. 445800 

R ( END  WALL  2 

)  = 

.  445800 

R  ( FLOUR 

)  = 

. 207900 

CONDITION  NUMBER  = 


2.  44? 


INTERNAL  WALL  SURFACE  TEMPERATURES  <C) 


1 ( ROOF 

T (SIDE  WALL  1 
T  (SIDE  WALL  2 
T (END  WALL  1 
T  (END  WALL  2 
T  (FLOOR 


4.55 

3.  76 
3.76 
3.76 
3.  76 


HEAT  LOSS  THROUGH  THE  WALLS  (W) 

Q  (ROOF  )=  314.91 

Q(SIDE  WALL  1  )=  157.46 

QtSIDE  WALL  2  )=  157.46 

Q (END  WALL  1  )=  97.91 

Q  (END  WALL  2  )=  97.91 

Q (FLOOR  )=  204.41 

RADIANT  HEAT  TRANSFER  BETWEEN  WALLS 


uRAD ( ROOF 

.-SIDE  WALL  1 

)  = 

6 . 80 

QRAD  (RiJOF 

.-SIDE  WALL  2 

)  = 

6.  80 

OR AD (ROOF 

.-END  WALL  1 

)  = 

3.96 

QRAD ( ROOF 

?  END  WALL  2 

)  = 

3 .  96 

QRAD ( ROOF 

-  FLOOR 

)  = 

50 .  96 

QRAD (SIDE 

WALL  1 

> SIDE  WALL  2 

)  = 

.  00 

QRAD (SIDE 

WALL  1 

.-END  WALL  1 

)  = 

.  01 

QRAD (SIDE 

WALL  1 

.-END  WALL  2 

)  = 

.  01 

QRAD (SIDE 

WALL  1 

-  FLOOR 

)  = 

10.  44 

QRAD (SIDE 

WALL  2 

-  END  WALL  1 

)  = 

.  0 1 

QRAD (SIDE 

WALL  2 

.-END  WALL  2 

)  = 

.  0 1 

QRAD (SIDE 

WALL  2 

-  FLO  OF; 

)  = 

10.44 

QRAD (END  WALL  1 

.-END  WALL  2 

)  = 

.  00 

QRAD (END  WALL  1 

-  FLOOR 

)  = 

6 .  02 

QRAD (END  WALL  2 

-  FLOOR 

)  = 

6.02 

TOTAL  HEAT  LOSS  BY  CONVECTION  AND  RADIATION  IS 


1030.05  (W) 


MEAN  THERMAL  INSULATION^ 
■♦'STOP* 


.053395  (K/W) 
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two  APC's  in  various  insulating  configurations.  The  measured  heat  loss  from  an 
uninsulated  APC  was  in  close  agreement  to  the  heat  loss  predicted  by  the 
theoretical  analysis.  However,  the  measured  heat  loss  from  an  APC  insulated 
with  a  covering  developed  by  the  Defence  Research  Establishment  Suffield  was 
found  to  be  significantly  greater  than  that  predicted  by  the  theoretical 
analysis.  Theories  are  advanced  on  possible  sources  of  error  and  discrepancy 
between  the  model  predictions  and  the  measured  values  of  temperature  and  heat 
flow.  Heat  loss  due  to  ventilation  is  also  examined  briefly. 
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